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Abstract. The Newton polytope Pf of a polynomial / is well known to have a strong impact on its zeros, 
as in the Kouchnirenko-Bernstein theorem on the number of simultaneous zeros of m polynomials with given 
Newton polytopes. In this article, we show that Pt also has a strong impact on the distribution of zeros of 
one or several polynomials. We equip the space of (holomorphic) polynomials of degree < N in m complex 
variables with its usual SU(m + l)-invariant Gaussian measure and then consider the conditional measures 
7|jvp induced on the subspace of polynomials whose Newton polytope Pf C NP. When P = S, the unit 
simplex, then it is obvious and well-known that the expected distribution of zeros Zt , (regarded as a 
current) of k polynomials /i, • ■ • ,/fe of degree N is uniform relative to the Fubini-Study form. Our main 
results concern the conditional expectation E\ NP (Z f 1 ,...,/ fc ) of zeros of k polynomials with Newton polytope 
NP C NpTi (where p = degP); these results are asymptotic as the scaling factor N — > oo. We show that 
E|jvp(^/ 1 ,...,/ fc ) is asymptotically uniform on the inverse image Ap of the open scaled polytope ip° via 
the moment map fj, : CP m — ► £ for projective space. However, the zeros have an exotic distribution outside 
of Ap and when k = m (the case of the Kouchnirenko-Bernstein theorem) the percentage of zeros outside 
Ap approaches as N — > oo. 
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Introduction 

It is well known that the Newton polytope P/ of a holomorphic polynomial f(zx, . . . , z m ) of degree p has a 
crucial influence on its value distribution and in particular on its zero set. Even the number of simultaneous 
zeros in C* m := (C \ {0}) m of m generic polynomials fx, . . . , f m depends on their Newton polytopes Pf j 



Be, Kol 



Ko2| . Our purpose in this paper is to demonstrate that the Newton polytope of a polynomial / 



also has a crucial influence on its mass density \f(z)\ 2 dV and on the spatial distribution of zeros {/ = 0}. 
We will show that there is a classically allowed region region where the mass almost surely concentrates 
and a classically forbidden region where it almost surely is exponentially decaying. The classically allowed 
region is the inverse image p^i^P) of the (scaled) polytope ip under the standard moment map /i£ of 
CP m . The simultaneous zeros of m generic polynomials fx, . . . ,f m in C* m almost surely concentrates (in the 
limit of high degrees) in the classically allowed region, giving a kind of quantitative version of the Bcrnstein- 
Kouchnirenko theorem. The Newton polytope has an equally strong (though different) impact for k < m 
polynomials fx, ■ ■ ■ , fk- The image of the zero set of fx , . . . , fk under the moment map is (up to a logarithmic 
re-parametrization) known as an amoeba in the sense of GKZ, Mi2fl . Results on the expected distribution 



of amoebas can be obtained from our results on the expected zero current; for example, as a consequence of 
Corollary^, for polytopes ^P with vertices in the interior of the standard simplex E C M." 1 , there is also a 
very forbidden region which the amoeba almost surely avoids. These patterns in zeros are statistical — they 
hold for random polynomials with prescribed Newton polytope — and are asymptotic as the degree of the 
polynomials tends to infinity. 

To state our problems and results precisely, let us recall some definitions. Let 

f(zi,...,z m ) = ^2 c a x a (zx,---,z m ), Xa(z) = z" 1 • • -z^ m 

aGN m :|a|<p 

be a polynomial of degree p in m complex variables. By the support of / we mean the set 

S f = {aeZ m :c a ^ 0}, (1) 
and by its Newton polytope Pf we mean the integral polytope 

P f := S f = the convex hull in E m of S f . (2) 



Throughout this paper we will make the simplifying assumption that P is a Delzant polytope [Dc, Gul] 



Non-Delzant polytopes are equally important in the study of polynomials, and we plan to extend our results 



to more general convex integral polytopes in the sequel | SZ2 . We also plan to study a parallel tunneling 
phenomenon for critical points. 

Our aim is to study the statistical patterns in polynomials with a fixed Newton polytope. We therefore 
define a natural Gaussian measure d^p on the space of polynomials with support contained in a fixed Newton 
polytope P. Since our purpose is to compare zero sets and masses as the polytope P varies, we view the 
polytope as placing a condition on the Gaussian ensemble of all polynomials of degree p equal to that of 
P, i.e. we define the measures d"/p as conditional probability measures of one fixed Gaussian measure. To 
define this ensemble, we homogenize the polynomial / to obtain a homogeneous polynomial F of degree p in 
m + 1 complex variables. Such polynomials may be identified with holomorphic sections F £ P°(CP m , 0(p)) 
of the p th power of the hyperplane section bundle 0cp>™(l). We denote by 

H (CP m ,O{p),P) = {Fe H (CF m ,O(p)) : P f C P} (3) 

the space of homogeneous polynomials F of degree d whose associated inhomogeneous form f(zx, ■ ■ ■ , z m ) = 
P(l, zx, ■ ■ ■ , z m ) has Newton polytope Pf contained in P. 

Let us recall that H (C¥ m ,O(p)) carries the standard SU(m + l)-invariant inner product 



(Fx,F 2 ) = [ FxF 2 da, 

J S 2m+l 



where da is Haar measure on the (2m+l)-sphere S 2m+1 . (See ( p8| ) for an alternate description in terms of the 
Fubini-Study form on CP™ and the Fubini-Study metric on 0(1).) An orthonormal basis of P°(CP m , 0(p)) 
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is given by {||xa|| 1 Xa}^ <p , where |j • |j denotes the norm in P°(CP m , 0{p)). (Note that ||x Q || depends on 
p.) The corresponding SU(m + l)-invariant Gaussian measure j p is defined by 

n P \\Xa\\ 
\a\<p 

where k p = #{a : \a\ < p} = ( m ^ p ). Thus, the coefficients A a are independent complex Gaussian random 
variables with mean zero and variance one. We then endow the space P°(CP m , 0(p), P) with the associated 
conditional probability measure 7 p |p: 

dl P \p{s) = ^p-e _|A|2 dA, s = ^2 Aq TTH[ ' ( 5 ) 

where the coefficients X a are again independent complex Gaussian random variables with mean zero and 
variance one. (By a slight abuse of notation, we let J2atP denote the sum over the lattice points a £ PnZ m ; 
#P denotes the cardinality of PnZ m .) We observe that, as a subspace of PT° (CP™, 0(p)), P°(CP m , 0(p),P) 
inherits the inner product (si,s 2 ) and that jip is the induced Gaussian measure. Probabilities relative to 
7|p can be considered as conditional probabilities; i.e. for any event E, 

Prob 7 {/ S E\P f =P} = Prob 7 | p (P). 

We now wish to consider the expected distribution of mass and zeros of polynomials with fixed New- 
ton polytope. ft turns out to involve the moment map /is of projective space. We use a notation for 
this moment map which reflects its dependence on the standard unit simplex S in ]R™ 1 with vertices at 
(0, . . . , 0), (f , 0, . . . , 0), (0, 1, . . . , 0), . . . , (0, . . . , 0, I ). The 'projective moment map' is then given by 

„ (-a ( N 2 l^l 2 

\i+\\z\r--'i+\\zf 

Via this moment map we define: 

Definition: Let P C be an integral polytope. The classically allowed region for polynomials in 
H (CP m ,O(p),P) is the set 

A P := fe 1 (-P°) C C" n 



(where P° denotes the interior of P), and the classically forbidden region is its complement C* m \ Ap. 

Our first result concerns the simultaneous zero set of m independent polynomials in m variables. Bezout's 
theorem tells us that m generic homogeneous polynomials Fi , . . . , F m of degree p have exactly p m simulta- 
neous zeros; these zeros all lie in C* m , for generic Fj. In fact, one immediately sees (by uniqueness of Haar 
probability measure) that the expected distribution of zeros is uniform with respect to the Fubini-Study 
projective volume form, when the ensemb le is given the SU( m + l)-invariant measure d^ p . 



According to Kouchnirenko's theorem [ Kol , Ko2 , At , Au |, the number of joint zeros in C* m of m generic 



polynomials {fi, . . . , f m } with given Newton polytope P equals m!Vol(P), where Vbl(P) is the Euclidean 
volume of P. For example, if P = pS, where E is the standard unit simplex in K. m , then Vol(pE) = 
p" l Vol(E) = E-j, and we get Bezout's theorem. (More generally, the fj may have different Newton polytopes, 
in which case the number of zeros is given by the Bernstein-Kouchnirenko formula as a 'mixed volume.') 

Now consider m independent random polynomials with Newton polytope P, using the conditional prob- 
ability (||). We let E\p(Zf u ...j m ) denote the expected density of their simultaneous zeros. It is the measure 
on C* m given by 

E|p(%,...,/J(C0 = / d lp]P {h) - J d JplP (f m ) [#{z £ U : h(z) = ■■ ■ = f m (z) = 0}] , (6) 

for U C C* m , where the integrals are over P°(CP m , Ok>), P). In fact, E| P (Z /l! ... Jm ) is an absolutely 
continuous measure given by a C°° density (see Prop. |4.2| ). Our first result is the surprising fact that as 
the polytope P expands, these zeros are expected to concentrate in the classically allowed region and are 
(asymptotically) uniform there: 
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1 f < 0n A P 

in the distribution sense; i.e., for any open U C C* m , we have 
1 



Theorem 1. Suppose that P is a Delzant polytope in K m . Then 

uj m or 

FS 

on C* m \A P 
e have 

-— -V lNP (#{z e U : h{z) = ••• = f m (z) = 0}) -» m!Vol CPm ([/n A P ) . 

In fact, our methods imply that the convergence of the zero current on the classically allowed region is 
exponentially fast in the sense that 

E|ivp(%,.,/J = (WX"s + O (e~ XN ) on A P , 

for some positive continuous function A on Ap. 

Theorem |l] is a special case of a general result stated below (Theorem |4|) on the expected distribution 
of simultaneous zeros of k random polynomials with given Newton polytopes P\, . . . , Pk- Although we only 
discuss expected behavior of zeros here, it is possible to show (using methods similar to those of [3Z1]) 
that distribution of zeros is self-averaging: i.e., almost all polynomials exhibit the expected behavior in an 
asymptotic sense. For the sake of brevity, we concentrate here on only the essentially new aspects of our 
problem, namely the dependence of the expected behavior on P. 

The following illustration^ shows the classically allowed region (shaded) and the classically forbidden 
region (unshaded) when P is the unit square in K 2 (and p — 2): 




Figure 1. The classically allowed region for P = [0, 1] x [0, 1] 
The terminology is taken from the semiclassical analysis of ground states o f Sc hroedinger operators Hh 



—h 2 A + V on K". The well-known 'Agmon estimates' of ground states (cf. |AgJ) show that £ 2 -normalized 
ground states or low-lying eigenfunctions H^tp — Etp of H are concentrated as h — > in the classically 
allowed region Ce = {x G K™ : V{x) < E} and have exponential decay |</?(x)| 2 = 0{e~ dtyX,CB ^ h ) as H — > in 
the complement. Here, oI(x,Ce) is the 'Agmon distance' to the classically allowed region, i.e. the distance 
in the metric y/V(x) — E dx. In our setting, the Hamiltonian is d*d on £ 2 -sections of powers 0(Np) of the 
hyperplane section bundle, H = 1/N, and the ground states are the holomorphic sections £T (CP m , 0(p), P). 



1 We would like to thank Aaron Carass for providing the illustrations in Figures |j) Q- [| an d 0- 
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Replacing |<^(a;)| 2 is the expected 'mass density' of the polynomials in this subspace. By definition, the mass 
density of a polynomial / at a point z G C* m is the Fubini-Study norm 

i/(*)i» := ( i+SV 2 (5=deg/) - (7) 

Out next result describes the asymptotics of the expected mass density of £ 2 -normalized polynomials 
with Newton polytope NP, i.e. the expected density Ev NP (\f(z)\l s ) with respect to the Haar probability 
measure, denoted by vnp, on the unit sphere in H°(Cf m , O(Np), NP). One easily shows (see §||) that 

HWI/WIL) = ^Wp-f^p (I/WIL) = #^pfW^) - ( 8 ) 

where H\np is the Szego kernel of H Q ( < CF m ,0(Np), NP). Here, the term 'Szego kernel' of a space S of 
holomorphic sections of a line bundle refers to the kernel for the orthogonal projection to iS from the space 
C 2 (dD) of square-integrable functions on the boundary of the associated unit disk bundle D; i.e., it is of 
the form H(x,y) — J^ . Sj(x)sj(y), where {sj} is an orthonormal basis of S. Precise definitions and the 



relationship to the classical Szego kernel for a strictly pseudoconvex domain are described in §1.2. Note that 

#{NP) = dim H (CF m ,O(Np),NP) = Yol{P)N m + . . . , 

which is a polynomial given by the Riemann-Roch formula (see §|3|). 

We shall show that the expected mass density E Vf/p (\f{z)\^, ) is asymptotically uniform with respect 
to Fubini-Study measure in the classically allowed region (as if there were no constraint at all); while in 
the forbidden region the mass decays exponentially. Thus in the semiclassical limit N — > oo, all the mass 
concentrates in the classically allowed region. To describe the behavior in the forbidden region, we subdivide 
it into regions corresponding to the faces of the polytope ip that lie in the interior S° of the simplex. To 

describe these regions, we recall that the fan associated to ^P (or P) is the collection of polyhedral cones 
C F C M m normal to the faces F of ±P (see § |l.l.l| ). For each face FcS°,we let 

TZ F = {e T -z-.t^Cf, /i E (z)ef}. (9) 



Note that if F is the open face ^P°, then Hp is the classically allowed region Ap. We show in § 3.2.1 that 
the classically forbidden region decomposes into the (disjoint) forbidden subregions 

jft F : F is a face of ±P, F C £°, F ^ ±P°j . 

We call IZf the flow-out of the face F. Boundary points of the TZf are called transition points. Points of 
dAp are always transition points, but there may be others. 

For the case where P is the unit square, then there are two interior faces whose flow-outs are the connected 
components of the classically forbidden region (see Figure [l] above). In this case, the set of transition points 
equals dAp. We shall give an example with an interior vertex, where the forbidden region is connected but 
decomposes into flow-outs of 3 faces (see Figure |?] in §5.2.1); in this case there are transition points not in 
dAp. 

Theorem 2. Suppose that P is a Delzant polytope in M. m . Then the expected mass density of random C 2 
normalized polynomials with Newton polytope NP has C°° asymptotic expansions of the form: 



%\Npt 3) = c ° + ^ N ~' + C ? N ~ 2 + ■ ■■ > for zeA P = ^ x (iP°) 

N -s/2 e -Nb{z) [ c F (z) + c F^ N -l + cF( z )N- 2 + ■••], forzell F 



where Cq and b\n° F are positive C°° functions on 1Z F , and s = codimP (for each face F C S° ). Furthermore, 
the remainder estimates in the expansions are uniform on compact subsets of the open regions 1Z° F and of 
A P . 

We shall use to restate Theorem || in terms of asymptotic formulae for the conditional Szego kernel on 
the diagonal for which we give more precise asymptotics in Proposition 3.1. In [ SZ2 we will extend the result 
to the non-Delzant case as well. Heuristically, the mass concentration in Ap can be understood as follows: 
the mass of a monomial m a concentrates (exponentially) on the torus ^^(a). The constraint a £ NP thus 
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concentrates all the mass of polynomials with fixed Newton polytope NP in Ap, with exponentially small 
errors, and the mass there is uniformly distributed with small errors. 
In particular, 



, for z e C* m \ Ap 

as illustrated by the following graphs, plotted using Maple, for the case where P is the unit square. (Recall 
Figure [l] for the depiction of Ap for this case.) 




Figure 2. %E UNP (\f(z)\ 2 FS ) for P = [0, 1] x [0, 1] 



The decay function b in Theorem || is analogous to the Agmon action to the allowed region. We shall 
show in Pro posi tion [3~l| that b is C 1 (but not C 2 ) on all of C* m , and we derive several formulas for b(z) 
((|l0|), (11), (|116| )). To describe our first formula, we associate with each point z £ C* m \ Ap a unique point 
£ G dAp of the form £ = e T//2 • z, where — r is in the (real) normal cone to the convex set P at the point 



p/^s(£) € BP. Here r 



(n^i, • 



denotes the 



action on 



1.1.1 



(The normal cone to a point 



is the element of the fan associated to the face of P containing q.) We write 



q € dP, described in 
t z = r, q(z) = these vectors are then given by the conditions: 

• PMs(e r * /2 • z) = q(z) e dP; 

• — t z is in the normal cone to P at q(z). 



The existence and uniqueness of q(z) and r z are stated in Lemma 3.8 
Our first formula for 6 is : 

i + IMI 2 



b(z) = (q(z),T x ) +p\og 



1 



oT z /2 



4 



We will obtain ( |io|) from the integral formula, 



b(z) 



q{e al2 ■ z) -pfi s (e 



a/2 



for zeC* m \ A! 



da 



(10) 



(11) 



(where Jq Z denotes the integral over any path in R m from to t z ), which could be interpreted as an action, 
thereby bringing the results closer to classical Agmon estimates. 

We also show that ^ logE|jvp (|/(^)| 2 S ) — > —b(z) uniformly on compact subsets of C*™ 1 (Proposition 3.2). 



RANDOM POLYNOMIALS WITH PRESCRIBED NEWTON POLYTOPE 



7 



Our final results concern the zero set of one or more random polynomials. We shall consider not only the 
expected volume density of the zero set, but also the more precise description involving the expected zero 
current. The zero current Zf lt ...j k of k polynomials (or more generally, k holomorphic functions) /i, . . . , fk 
is the current of integration over the zero set 

\Zh,.,f>\ i z e C* m : h(z) = ■■■ = f m {z) = 0} . 

This current is given by 



(% fk ,cp):= ip, for a test form tp e T> m ~ k > m - k (€* m ) , 

J\z fl ,..., fk \ 

where integration is over the set of regular points of \Zf lt ...j k \ (provided that it is of codimension k and 
without multiplicity, which is almost surely the case). 

We denote by E|p(Zy) = E 7p |p(Zj) the conditional expectation of the zero current of a random polynomial 
/ e H (CP m ,O(p),P) with Newton polytope P. In fact, E\ P (Z f ) is actually a smooth (1, l)-form on C* m 
(Proposition |4.l| ) . 

Let us recall what happens when P — pE (cf. JSZ1 1 ) . By the uniqueness of Haar measure, the expected 



zero current ~Ei{Zf) taken over all polynomials of degree p is given by pu> FS , where aj FS = ^dd log || z\\ 2 is the 
Fubini-Study Kahler form on CP m . Thus the expected distribution of zeros, as well as the tangent to the 
zero varieties, is uniform over CP m . We now consider how the expectation changes if we add the condition 
that P f = P. 

Theorem 3. Let P be a Delzant polytope. Then there exists a closed semipositive (1, l)-form ipp on C* m 
with piecewise C°° coefficients such that: 

i) N- l B\ NP {Z } )^i> P mC\ oc (C* m ). 

ii) ipp = pu> FS on the classically allowed region ^^(iP ). 

iii) On each region 1Z° F , the (1, l)-form ipp is C°° and has constant rank equal to dimP; in particular, if 
»6S° is a vertex of ^P ; then ipp\n° = 0. 

As a corollary, we obtain some statistical results on the so-called 'tentacles' of amoebas in dimension 2 



(see p.3|). Roughly speaking, the (compact) amoeba of a polynomial f(zi, z?) is the image of the Riemann 
surface Z/ under the moment map /i£ on (C*) 2 , and the tentacles are the ends of the amoeba. Certain 
tentacles must end at vertices of the triangle E while others are 'free' to end anywhere along the boundary 



of E. In Corollary 4.7, we will prove that (in the limit N — > oo) almost all of the free tentacles of typical 
amoebas tend to end in the classically allowed portion of <9E. 

We call the form ipp in Theorem || the limit expected zero current. By C\ oc convergence in (i), we mean 
£ 1 1 oc convergence of the coefficients. (Recall that E|jvp(^/i,. ..,f k ) is a (k, fc)-form with smooth coefficients.) If 
we write ipp = \/— 1 ^2 i^jki^dzj Adzfc, then (ipjk(z)) is a semi-positive Hermitian matrix, for non-transition 
points z. By the rank of ipp at z, we mean the rank of the matrix (ipjk(z)). Note that if TZf and TZpi are 
adjoining regions (i.e., have a common codimension 1 interface), then F and F' are of different dimensions, 
so ipp must be discontinuous along the interface. We remark that as an element of I?' 1 ' 1 (C* m ), the current 
ipp is closed and positive. 

The form ipp not only encodes the expected (normalized) density of the zero set, but also the expected 
density of tangent directions to the zero set. In the course of the proof of Theorem ||, we will show that in 
the forbidden region, the limit tangent directions are restricted. In particular, as the polytope expands, the 
tangent spaces to typical zero sets through a point z° € 1Z° F approach spaces containing the tangent space 
(at z°) of the 'normal flow' {e T+l0 ■ z° : t,9 E Tp C W n }. Precise formulations of this fact are given in 



Theorem 4.5 and (162). Thus, while the expected distribution of zero densities is absolutely continuous, the 
expected distribution of zero tangents is singular. 

We can also consider k independently chosen random polynomials 

/^^(cr.ofe),^), i<j<k 
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(k < to), and we let E|p l! ....p k {Z f 1 ....j k ) denote the expected zero current with respect to the probabil- 
ity measure 7 Pl |p 1 x ••• x J Pk \p k on the product space. If Pi = • • • = Pk = P, then we also write 
E \p,...,p( z h.---Jk) = E \p( z h,---Jk)- 

Theorem 4. Let P\,...,Pk be Delzant polytopes. Then 

N- k V ]NPl _ NPk {Z fu ...j k )^ip Pl /\---A^ Pk m 4(C™), as N ^ ^ . 



We see from Theorem |^(iii) that the zero set \Zf\ of a polynomial with polytope NP almost surely creeps 
into the classically forbidden region //^ (S \ |P) in the semiclassical limit N — > oo. Indeed the expected 
volume of the zero set, or more generally the simultaneous zero set of k polynomials, has the following exotic 
distribution law: 

Corollary 5. Let P%, . . . , P/. be Delzant polytopes. Then for any open set U CC C* m , 

^^,...,^01(1%,...^ | D 17) - ^^yy /^fiA-A 4>P k A ■ 

By the volume of f h \ (1 U, we mean the (2m — 2/c)-dimensional volume in CP m , which is given by 

Vol(|%,..., /fc | n U)= f j^^ m ~ k = (%,...,/*, T^w^ k ) ■ ( 12 ) 

J \z fl fk \nu v ' 

Corollary |^ follows immediately from Theorem |i] and ( |l2| ) . 

In addition, the following consequence of CorollarjTjHI and part (iii) of Theorem || says that there are 
subtler 'very forbidden regions' in the case where the polytope has faces of codimension > 2 in the interior 
of S, namely the regions comprising the normal flow-out of these faces: 

Corollary 6. The expected zero current N~ k ~E\ NP (Zf 1 : ...j k ) tends to at all points of each forbidden 
subregion IZf with dimF < k. In particular, N^ 1 'E\ i ^ P {Z f) tends to zero at all points of the flow-out 1Z V of 
a vertex v G ~P n S° . 

We now say a few words about the proofs. The key result is Theorem p| o n the mass of polynomials 
with Newton polytope P, which we formulate more precisely in Proposition 3.1 using the conditional Szego 
kernel. As will be shown in we can derive the expected distribution of zeros from our asymptotic formula 
for the conditional Szego kern el Tl\ PN (z, z). 

As in the proof in [At, Au| of Kouchnirenko's theorem, a natural approach to our results on polynomials 
with prescribed Newton polytope is through the Kahler toric variety (M p ,uj p ) associated to P. We recall 
that M P carries a natural line bundle L P with c\{L P ) — ui P . The main connection is that there exists a 
natural identification 

H (CP m ,O(Np),P) ~ H°(Mp,L$) (13) 

between polynomials with prescribed Newton polytope and holomorphic sections of powers L p of L P . The 
torus T m = {(e l(pi , . . . , e^™)} acts on H°(M P , L P ) with character 

XNp(e iip ) = V , e* = (e^\. . .,e^) . (14) 



aeNP 



We then have the simple expression 



U\np{z,z)= Tt% {t-z,z) X Np{t)dt. (15) 



To obtain asymptotics as N — > oo, we derive an oscillatory integral formula. First we combine ( |15| ) with the 
formula 

XNp(e lv )= f n^(e^ ■w,w)dVo\ Mp H, (16) 

JMp 

where Ii Mp is the Szego kernel of M P , i.e. the orthogonal projection onto H°(M P , L p ) with respect to the 
volume form determined by uj p . In order to obtain an exponential decay rate e~ Nb ^ of Theorem 0, we will 
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need a very precise understanding of H^ Ip 



i x . This we obtain in §|| by making a special construction of the 
Szego kernel of a toric variety which is perhaps of independent interest. Our formula for II p involves the 
application of a certain Toeplitz Fourier multiplier to t he c lassical Szego kernel for the sphere, pulled back by 
a monomial embedding. Aft er showing in Proposition 2.2 that this Toeplitz- Fourier multiplier is a Toeplitz 
operator in the sense of BG], we obtain an oscillatory integral formula for H\np (see (|69| ) and (ff2|)) in terms 
of an integral over T m x Mp with complex phase of positive type. 

After that, the proof of Theorem || follows from the method of stationary phase for complex oscillatory 
integrals. For the case where z is in the classically allowed region, we easily find that the critical manifold 
is given by tp — and /ip(w) = pfis(z), where [ip is the moment map of Mp. Since the phase vanishes 
and has nondegenerate normal Hessian along the critical manifold, we immediately obtain an asymptotic 
expansion. The case where z is in the classically forbidden region is more subtle. Since pfi^(z) lies outside 
of P for this case, the phase has no critical points. To complete the analysis, we must deform the contour to 
pick up critical points. The main part of the proof of Theorem |^ involves finding a contour for which there 
is a nondegenerate critical manifold where the phase has maximum real part. In particular, we consider the 
complexification C* m of T m and deform T m to a contour of the form (log . . . , log \ C, m \) = t € K m . We 
then show that the phase takes its maximal real part on a critical manifold if and only if r is the unique 
vector t z used in formula (^) for the decay rate b(z); indeed the maximal real part of the phase is — b{z). 
Finally we prove that nondegeneracy of the normal Hessian holds if (and only if) z is not a transition point, 
so that we obtain the asymptotic expansion of Theorem ^. 

We would like to remark here that other formulae for the lattice sums xp(t) are well known from the 
work of Khovanskii-PukhlikoyJ KP | , Brion-Vcrgne JBV| and others. An alternative approach is to replace 
( |l6| ) with the formula from [BV, Theorem 3.12]: 



Xp(e ilp ) = Todd(<9/<%) 



<p(h) 



(17) 



h=0 



where P(h) = {x : (uj, x) + dj + hj > 0, 



1 < j < n }, and Todd(d/dh) is a certain infinite order differential 



operator known as a Todd operator. Upon dilating the polytope, one obtains 



XNp{e iip ) = N m Todd(Ndfdh) 



e iN(x, V ) dx 



lp(h) 



(18) 



h=Q 



The asymptotics of mass and zeros in the classically allowed region follow very easily from this expression 
for xnp, and with further effort one can prove the result in the forbidden region. Both approaches seem to 
us quite interesting and to have their own advantages. In the sequel [3Z2|, we will use the approach based 
on (fl8|), which seems more efficient when the toric varieties have singularities. 

We end the introduction with some comments on the relation of our results to other work on polynomials 
with a fixed Newton polytope. Such polynomials are called sparse in the literature, and methods of algebraic 
(including toric) geometry have recently been applied to the computational problem of locating zeros of 
systems of such sparse polynomials (e.g., see [HS, [MR| , [5t| , [Ro| > |Ve[). Our results give information on the ex- 
pected location of zeros when the polynomials are given the conditional measure J\np on H°(CP m , O(N), P). 
To our knowledge, the asymptotic concentration of zeros of sparse systems in the classically allowed region 
has not been observed before. It obviously pays most to search for zeros of systems of m polynomials in m 
variables in the allowed region. 

It should be noted that this asymptotic pattern reflects the choice of the conditional measures j\np on 
H°(CP m , O(N), P). In previous work [SZ1|, we obtained the expected and almost everywhere distribution of 
zeros of several holomorphic sections of positive line bundles L — > M over general Kahler manifolds. Those 
results specialize to our current setting where M = Alp, L = Lp and also give results on the distribution of 
zeros of polynomials in i?°(CP m , O(N), P). However, in that paper we defined different Gaussian probability 
measures ^f N F on H°(Mp, Lp), namely those induced by the C 2 inner product induced by a hermitian metric 
on Lp and its curvature form up. In Proposition l.£ of [ 5Z1| , we showed that the expected distribution of 
zeros relative to the r y^ p satisfies ir^E M P (Zf lt ...j k ) — up + O(j^). The measures have more recently 



been studied by MR|, who also obtain (among other things) the formula for the expected distribution. Our 
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present measures 7^^ are quite singular relative to 7|jvp in the limit as N — > oo. Our point of view is that 
the conditional measures 7|jvp on i?°(CP m , 0(N), P) are the natural ones when comparing polynomials in 
H° (CP m , O (N) , P) as P varies. The toric variety Mp is a technical device for studying II|jvp but does not 
figure into the statements of our main results. The difference in the two measures is that in the conditional 
case, we fix the C 2 norms of the monomials z a in advance (as their CP™ norms), while in the 7^f p case the 
norms vary as P and Mp vary. The latter would seem to create complicated biases towards some monomials 
and away from others as P varies and make it difficult to understand what such a comparison is measuring. 



1. Background on toric varieties and moment polytopes 

1.1. Newton polytopes and toric varieties. The space H°(CP m , G(j>), P) of polynomials with Newton 
poly tope P may be identified with the space H°(Mp, Lp) of holomorphic sections of a line bundle Lp — > Mp 
over a toric variety Mp. Under our running assumption that P is Dclzant (see below), the variety Mp is 
smooth. Recall that a toric variety is a complex algebraic variety M containing the complex torus 

C* m = (C\{0» x ••• x (C\{0» 

as a Zariski-dense open set such that the group action of C* m on itself extends to a C* m action on M. In 
the Delzant case, Mp can be given the structure of a symplectic manifold such that the restriction of the 
action to the underlying real torus 

T m = {(Ci, ■ • ■ Cm) G C* m : |0| = 1, 1 < j < m} 



is a Hamiltonian action (see §1.3). 



We now review two (well-known) constructions of Mp, by fans, resp. monomial embeddings, to establish 



notation and clarify the properties we will be needing. For further details, see [Da, Fu, Oda|. We point 



out here that there is also a geometric approach due to Delzant [De] through symplectic reduction; see 



[De, Au, Gul] 



1.1.1. Fans. By a convex integral polytope, we mean the convex hull in K n of a finite set in the lattice 1 m . 
A convex integral polytope P with n facets (i.e. codimension-one faces) can be defined by linear equations 

£i(x) := (x, m) + en > 0, (i = 1, . . . , n), 

where u, € Z m is the primitive interior normal to the i-th facet. The polytope P is called Delzant if each 
vertex is the intersection of exactly m facets whose primitive normal vectors generate the lattice Z m . 
For each point x G P, we consider the normal cone to P at x, 

C x := {u £ W n : (u,x) = sup(w, y)} , 

which is a closed convex polyhedral cone. We decompose P into a finite union of faces, each face being an 
equivalence class under the equivalence relation x ~ y C x = C y . For each face F, we let Cf denote 

the normal cone of the points of F. Note that by our convention, the faces are disjoint sets. We 
shall use the term closed face to refer to the closure of a face of P. 

Each face of dimension r (0 < r < m) is an open polytope in an r-plane in R m ; i.e., the 0-dimensional 
faces are the vertices of P, the 1-dimensional faces are the edges with their end points removed, and so forth. 
The facets F^ and their normal cones are given by: 

Fi = {x e P : ii{x) = 0} , C Fi = {-t Ul : t > 0} . 

The m-dimensional face is the interior P° of the polytope with normal cone Cpo = {0}. 

A convex integral polytope P determines the fan Tp := {Cf ■ F is a face of P}. A fan T in R m is a 
collection of closed convex rational polyhedral cones such that a closed face of a cone in T is an element of T 
and the intersection of two cones in T is a closed face of each of them. An example of an integral polytope 
and its fan is given in the illustration below. 

There is a natural way to construct a toric variety from a fan by gluing together the affine varieties arising 



from the cones in the fan. As we will not use this construction here, we refer the reader to [Fu] for the 
details. We note that a fan whose union is all of W n defines a complete toric variety, which is projective if 
and only if T is the fan of a convex polytope P. If P is Delzant, then the toric variety of the fan Tp is a 
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Figure 3. A convex polytope and its fan 

smooth projective manifold. The polytope in Figure [| is Delzant if W3 = (0, 1), and the corresponding toric 
variety Mp is a Hirzebruch surface (see We shall study this example in §5.2.1. 



1.1.2. Monomial embeddings. An alternative definition of Mp is through T m -equi variant monomial embed- 
dings. We let P n Z = {a(l), a(2), . . . a(#P)}. If we fix c = (a, . . . , c # p) G (C*) #p , we get a map 

$p = [c Q( i)Xa(i), • ■ .,c a(#P) Xa(#p)] : C* m -» CP #p - 1 ; 

i.e., 

$p(*) = [ Cq(1) z«W..., Cq(#p) z«(# p )] . 

The closure of the image is a toric variety Mp C CP# P_1 . If P is Delzant, then $p is injective and we 
identify C* m with its image (the 'open orbit') in Mp. 
We refer to the resulting embedding 

$p : Mp <^> CP #P ^ (19) 

as a monomial embedding, and we define the line bundle 

hp := $£0(1). (20) 

Furthermore, H°(Mp, Lp) = $p* P°(CP# P_1 , 0(1)); i.e., the sections are spanned by the monomials \a- It 
follows that 

H°{Mp, L P ) ~ P°(CP"\ 0(p), P) . (21) 



We further have (see Fu | ) : 

H°{Mp, Lp) ~ $p*i/°(CP #p - 1 , 0(AT)) = i/ °(CP m , 0(Np),NP) . (22) 

Recall that P°(CP #P ~\ 0(1)) has as a basis the linear coordinate functions A, : C #p -> C, 1 < 7 < #P. 
The Fubini-study metric h FS on 0(1) is given by 

101 



|Ailps([C]) = g (Ce 

which has curvature form w FS = ^Sdlog ||C|| 2 - We endow Lp with the Hermitian metric hp := $p*/j FS of 
curvature cjp given on C* m by 

c P = $ c pX s = ^^log J] |c a | VP- (23) 

aeP 

Each monomial with a E P corresponds to a section of H°(Mp, L c p ) and vice versa. To explicitly 
define this correspondence, we make the identifications (recalling (|20|)): 

Xa(j) = c «tj)®P*^ e H °( M p> L< p)> 1 ^ 3 < #P ■ (24) 



12 



BERNARD SHIFFMAN AND STEVE ZELDITCH 



So far, we have not specified the constants c a . For studying our phenomena, the choice of constants 
defining the toric variety Mp is not important. However, when our polytope P is the full simplex pS, we 

1/2 

shall use the special choice c Q = ( p ) , where ( P J is t he multinomial coefficient defined in the next section. 
(We use c a = 1 in §||, but except for Proposition which is not used in our proof, any choice of c will 
work.) 

Remark: For a fixed polytope P, the symplectic manifolds {Mp,wp) are equivariantly symplectically equiv- 
alent, by Delzant's theorem [De, Gul, Au|. The Mp are also equivariantly holomorphically equivalent. 
How ever, in general the maps defining the equivalences are not simultaneously holomorphic and symplectic 



(see |Guf 



f .2. Szego kernels. We first describe a general Szego kernel for a positive Hcrmitian line bundle (L, h) on 
a compact complex manifold M. Let Sn C H°(M, L n ) denote a subspace of the (finite dimensional) vector 
space of holomorphic sections of the iV th power of L. We give Sn an inner product. An important example 
is the 'induced inner product': 

{si,s 2 )s N = (s 1 (z),s 2 (z)) dVol M (z) , s 1 ,s 2 eS N - (25) 
Jm x ' h N 

By the 'Szego projector' Hs N , we mean the orthogonal projection onto Sn- As in [Zc, 3Z1, BSZfl , we lift 
our Szego kernels on the associated principal S 1 bundle X — > M, defined via a Hermitian metric h as follows: 
Let L* — > M denote the dual line bundle to L with dual metric h*, and put X — {v e L* : \\v\\h* — 1}. We 
let r$x = ex (x € X) denote the S 1 action on X. We then identify sections sn of L N with equivariant 
functions §n on X by the rule 

§n(X) = (\® n ,s n (z)) , XeX z , (26) 

where \® N = A® • • • ® A; then sjv(^e^) = e lNe §n(x), and we regard elements of Sn as equivariant functions 
on X. 

If {-Sf} denotes an orthonormal basis of Sn, then the projector Hs N is given by the kernel 



U SN (x,y) =Y,S»(z)S»(y). (27) 
i=i 

(Note that the usual Szego kernel for the strictly pseudoconvex boundary X is given by X]jv=i ns N , where 
S N =H°(M,L N ).) 

We now describe three different sequences of Szego kernels on toric varieties, which play a crucial role in 
our main results: 

f .2.1. The projective Szego kernels. As a first example, we let L be the hyperplane section bundle 0(1) on 
CP m . Recall that the space H°(CF m , 0{p)) of holomorphic sections of L p = (D(p) consists of homogeneous 
polynomials 

F«0, ■ ■ ■ , Cm) - £ C X ( X (C A - Co A ° • • • Ct ) 

|A|=p 

in m + 1 variables. We can identify F with the (non-homogeneous) polynomial in m variables, 

f(zi, ...,z m ) = f(i,zi, ...,z m )= ^2 c <* za ( z ° = z T ■ ■■zZr), 

\a\<p 

where c a — C^p, c\ p = (p — \a\, a%, . . . , a m ). Inversely, a polynomial f(z) of degree < p can be identified 
with an element of #°(CP m , <D(p)). 

We give CP m the Fubini-Study Kahler form given in homogeneous coordinates (Co, ■ ■ ■ ,(m) by w FS = 
log ||C|| 2 , and we give 0(p) the Fubini-Study metric: 

\F(C)U = l^(C)l/IICII P , for F E H Q (CP"\0(p)) . 
Identifying F with the polynomial f(z) = F(l, z±, . . . , z m ), the Fubini-Study norm can be written 

\f(z)\ FS = \f(z)\/(l + \\z\\ 2 y/ 2 (zeC m ). 
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We equip the space S p = H (C¥ m ,O(p)) of all homogeneous polynomials of degree p with the inner product 

(/(*),$£)> 



(f,9) = [ (F>G) dVo\ Cpm = — f jf! 
J crm m \ J cm (1 



■«&(*), /, 5 eP u (CP"\0(p)) 



(28) 



+ lkll 2 ) p 

(We use here the Riemannian volume dVolcpm = ^r^^ ; note that the total volume of CP" 1 is ^ , using our 
conventions.) 

Under the T" 1 action, we have the weight space decomposition 

H (CP m ,O(p))= C Xq , 
\a\<p 

where we recall that Xa{z) = z" 1 ■ ■ ■ z," m . The monomials {xa} are orthogonal but not normalized. Any 
choice of norming constants {r a e C} will give a monomial basis {r a x a } for H (CP m ,O(p)). We shall 
choose r a = ,, ,, . where 



\\XaW = \J (Xa,Xc 



(p + m)\(l) 



[p- |a|)!ai! ■ ■ ■ a m \ 



is the Fubini-Study C 2 norm of Xa given by (p^). (See SZ1 , §4.2]; the extra factor ml in [3Z1 is due to the 
use of oj m instead of for the volume form.) This choice provides an orthonormal basis for H°(C¥ m , 0(p)) 
given by the monomials 



1 



\Xa\ 



' Xa 



(p + m)\ ( p 



p\ 



\a\ < p 



We can identify L* = 0Qpm(— 1) with C m+1 with the origin blown up, and the circle bundle X C L* is 
unit sphere S 2m+1 C C m+1 . We let x a : S 2m+1 -> C denote the equivariant lift of x Q € i/°(CP m , £>(p)). 
We note that 

Xq(x) = x qP , 

and hence the Szego kernel IIp P for the orthogonal projection is given by: 



i 



(p + m)\ 



Xa(^)Xa(2/) = 
|«|<p llXa " P ' |a|<p 

for x,y £ S 2m+1 . (The sum Hp 1 "" is the usual Szego kernel for the sphere.) 



(29) 



1.2.2. The conditional Szego kernels associated to a polytope P. In this case, the relevant space of polynomials 
is the subspace H°(C¥ m , 0(p), P) C H°(CF m , 0(p)) of polynomials with Newton polytope P. Here, we may 
choose any p > degP := max{|a| : a G P}, but we normally choose p — degP. In the conditional Szego 
kernel, we retain the Fubini-Study inner product, denoted (,)|p, on this subspace. Hence this example is 
very similar to the previous one. The main difference is that under the T m action, we have the weight space 
decomposition 

P°(CP™,O(p),P) = 0C Xq . 

Definition: The conditional Szego kernel Tl\p is the kernel for the orthogonal projection to P°(CP m , 0(p), P) 
with respect to the induced Fubini-Study inner product: 



R\p(z,y) = E I, 1 u 2 Xa(x)Xa{y) ■ 



(30) 



When defining the term 'random polynomial with fixed Newton polytope P', we wish to use an £ 2 -norm 
on monomials which is defined independently of P. This explains why the conditional Szego kernel is the 
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essential one in our problem. The conditional Szego kernel can be written explicitly as 



e ip(0-<p) ( p )z a w a 



where (z, 9), resp. (w, ip), are local coordinates for x, resp. y, in X; i.e., 



x = e ie (l + \\z\\ 2 )- 1 / 2 (l,z 1 ,...,z m ). 



1.2.3. The intrinsic Szego kernels of a toric variety. We now consider Szego kernels which are defined by 
intrinsic inner products associated to the toric variety Mp = Mp. (To simplify notation, we now drop the 
c.) We regard them as a technical device for obtaining properties of the conditional Szego kernels. 

We thus put iS/v = H°(Mp, Lp), equipped with the inner product (,)m p derived from a Hcrmitian 
inner product hp on Lp. As on all S 1 bundles associated to line bundles, the S 1 action on Xp gives the 
decompositions 

£ 2 (Xp) = ®^ 1 £%(X P ), feC% ^ f(e iS -x)=e im f(x). 
Restricting to the Hardy space, we obtain 

CO 

H 2 {X P ) = ($H 2 N {Xp), H 2 N = C 2 N n Ti 2 . 

N=l 



Recall that we have the canonical identifications TL 2 N (Xp) = H°(M P , Lp) given by identifying a section sjy 
with the equivariant function §n on Xp. We let (,)x P denote the corresponding inner product on Ti N (Xp), 
and we consider the orthogonal projectors Ii ! ^ p {x,y) relative to this inner product. 

We let Xa denote the equivariant lift of Xa to Xp. Since (,)x P is invariant under the T m action on Xp, 
the {Xa} are orthogonal: indeed, 

= (t*xZ,t%)x P = t a t^(x^xJ)x P , for all t e T m , 

and hence (x« , Xp)x P = unless a = [3. We then have: 

n%"(x,y)= u-^€(*)W¥), (32) 

aeNP WXaWxr 

where ||Xa||x P is the C 2 norm of Xa with respect to the inner product (,)x P - 

In particular, we note that My, = M p y, = CP™ 1 , and we have the circle bundles X p z — ► CP m , for p > 1. 
When p=l,X s = S 2m+1 while for p > 1, it is the lens space X pT , = 5 2m+1 /{e 27ri/p }. The latter statement 
follows from the fact that homogeneous polynomials of degree p are well defined on (and separate points of) 
the quotient by the cyclic group of p th roots of unity. 



1.2 .4. Powers of a line bundle. Our main results are asymptotics results in powers of a line bundle or 
(equivalently) dilates of the polytope P. We therefore take powers of the line bundles in the previous 
three examples obtaining the spaces H° (CP m , O (Np) , NP) and H°(M P , Lp) ~ H a (M NP ,L NP ). Hermitian 
metrics on the original bundles induce inner products on these spaces. 
To summarize, we have the following three sequences of Szego kernels: 

• n^ p ™ = the orthogonal projection to H a (CP m ,0(Np)); 

• Tl lNP = the orthogonal projection to H (CF m ,O(Np),NP); 

• n^ /p = the orthogonal projection to H°(M P ,L^). 

We note that for the case P = pS, the three sequences coincide. However, if P^pS, then the second and 
third sequences are quite different. 
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1.3. Moment maps and torus actions. The group C* m acts on Mp and the subgroup T m acts in a 
Hamiltonian fashion. Let us recall the formula for its moment map /ip : Mp — > R m , restricted to the open 
orbit C* m . This moment map is the composition 



//',. : Mp ^-- ^-i po 



where 



(j, ([zi,...,z #P }) = -^(\zi\ 2 ,...,\z #P \ 2 ), 



\4 

and A is the linear projection is given by the column vectors (a 1 , . . . , a #p ). Hence we have 

VP = V \l |2| Y 12 E M'lXal 2 " . (33) 

For any c, the image of Mp under \xp equals P. 

For a G P, we again lift Xq € H°(Mp, Lp) to an equivariant function % p on the circle bundle Xp — > Mp, 
and we write 

™a(j) := C «0)XaO) = Cj ° «-P (34) 

where ip : Xp — > S* 2d+1 is the lift of the embedding Mp °-> CP# P_1 . (Of course, m p depends on c, which 
we omit to simplify notation.) We also consider the monomials 

7Tl a . — CaXa 

so that m p is the equivariant lift of m a to Xp. In terms of local coordinates (z, 9) on 7r _1 (C* m ) C Xp, we 
have 

*) - — ; 1/2 • (35) 



(E, £ pM 2 M 2 ) 



Noting that 



we obtain the formula: 



#p 

£|tf£l 2 = £lCi°^l 2 = i> ( 36 ) 

Mp(*) = £ |m p (*)| 2 a . ( 3? ) 

(We write |m p (z)| = |m p (z, #)|, since the absolute value is independent of 6.) 

In the case where the polytope is the full simplex pS, we shall use the special choice 

i i 
p\ 2 ( p\ \ 2 n n_i 

' .. I 

so that 



MpS (z) := < s (z) = E r = 1 + f ■ |2 (H 2 , ■ ■ ■ , M 2 ) , (38) 

where the last equality follows by differentiating the identity (1 + J2xj) p — J2\ a \< p (a) 2 -"- Note that this 
choice gives us the scaling formula 

Furthermore, 

' a {z)l -[(m+p)l\ \\ Xa \\ ' (39) 
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where, by abuse of notation, we regard \\a\ and |m^ s | as functions on C* m , since they are invariant under 
the circle action. Recalling that M p s = 5 ,2m+1 /Z p , we have the more precise formula 



pi 



(m + p)\ 



\\Xa\\ 



where x' G M p s is the equivalence class of x G S + , and therefore by (|2Sf), 



£ mf (*')m?V) ■ (40) 

We now recall a fact on zeros of eigensections which will be useful later on. 
Proposition 1.1. Let a G P. Then for all z G M P , we have 

Xa(z) = Ltp(z) e|J{F:Fisa facet , a g F} , 
Proof. We first verify the implication (<=): Let z G fj,p (F), a £ F. The moment map (37) exhibits (J>p(z) 



as a convex combination of a G P. Points on F can however only be convex combinations of lattice points 
a G F. Hence the coefficients of a $ F must vanish. Therefore Xa( z ) = 0. 

Since the image under the moment map of the support of Div Xa is a union of facets of P, to prove the 
reverse implication, it suffices to show that for if F is a facet of P, then \ a ^ on /^p 1 (F) whenever a G F. 
But since /zp is surjective, we see from (^7|) that this is true whenever a is a vertex of F. (Otherwise, points 
of F near a would not be in the image.) Now suppose that a G F is arbitrary. We can write n^a — J^" 1 UjVj, 
where the Vj are the vertices of F, rij G N and Yl n j = no > 1. Then Xa° = Xvl ' ' ' XvZ ^ on Mp 1 ^)- C 

Remark: Recall that the polytope P is defined by the equations tj > 0, where £j(x) = (uj, x) + a,j, with Uj 
the primitive inward normal to the closed facet Fj — {£j = 0}. A more precise statement of Proposition |l.l| 
is 

Divxa = y^Jj(a)Dj , where D 3 = ^(Fj) , 

3 

for all a G Z m . In particular, letting a = 0, we see that Lp is the line bundle of the divisor Dp = ^ QjDj- 
We now review how the action of the real torus T m lifts from Mp to Xp and combines with the S 1 action 



to define a T m+1 action on Xp. We use the monomial embedding approach of § 1.1.2 : Recall that under the 
monomial embedding 

$ P : C™ - Mp - CP# p -\ z i * [c q(1) z q(1) , . . . ,c Q(#P) z Q (# p )] , 

the T m action on Afp C CP #P " 1 is given by 



e 



[Ci,...,C#p]= [e l<Q(1) ^Ci,---,e^# p )^C # p] . (41) 

The action ( ]4l| ) lifts to an action on Lp 1 : 

e^-C= (e^W'^Ci,--- ! e i<a(#p) ^ > C#F) . (42) 

Since the circle bundle Xp C S 2 ^ p ^ 1 is invariant under this action, ([42] ) also gives a lift of the action ( ^l| ) 
to X p . 

We also have the standard circle action on Xp: 

which commutes with the T m -action (f42|). Combining ( |42] ) and (|43|), we then obtain a T m+1 -action on Xp: 

(e ie ,e lVl ,...,e I¥,m )»C = e ie (e l¥, -C) • (44) 
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The torus action on Xp can be quantized to define an action of the torus as unitary operators on H. 2 (Xp) = 
©jv=o T^n(^p)- Specifically, we let Si, . . . E rn denote the differential operators on Xp generated by the T m 
action: 

(Z j S)(0 = - i -^S(e i f- 01^=0, SeC^(Xp). (45) 

Proposition 1.2. For 1 < j < m, 

(i) Z r .H%(Xp)^H%(X p ); 

(ii) The lifted monomials Xa £ TL 2 N (Xp) satisfy ^jXa — a jXa ( a € NP). 

Proof. Item (i) follows from the fact that the T m action is holomorphic and commutes with J^. For the 
case N = 1, (ii) follows immediately from ( |34| ) and (f42|). For a € NP, N > 1, we write a = /3 1 • • • (3 N with 
/3 fc e P, and the conclusion then follows from the first case and the product rule. □ 

Furthermore, we recall that 

— :?4(Xp) ^H 2 N (X P ) , ~—§ N = Ns N for sjv G TL 2 N (Xp) . (46) 



Remark: The vector fields Sj can be constructed geometrically as follows (see [Gul|): Let £j = -J^— (1 < j < 



m) denote the Hamiltonian vector fields generating the T m action on Mp. There is a natural contact 1-form 



a on Xp determined by the Hermitian connection; a key property of a is that da — ti*u (see [ |Zq , 3Z1 ). We 
use a to define the horizontal lifts of the Hamilton vector fields ^f. 

■K^ = i h <*($) = o. 

The vector fields Sj are then given by: 

(Here, ^* S K m is the element of the Lie algebra of T m which acts as £j on Mp.) 

2. SZEGO KERNELS ON TORIC VARIETIES 

As mentioned in the introduction, our analysis of the conditional Szego kernel is by means of the integral 
formula (|l6|) relating it to the intrinsic Szego kernels nj^ p of Mp. The purpose of this section is to give a 
special construction of the Szego kernels of a toric variety which is sufficiently precise for our applications. 



Our first step, Proposition |2.1|, is to give an exact formula for H^j p as the composition of a certain 



Toeplitz- Fourier multiplier denoted (VQ)^ 1 and the pull back of the Fubini-Study kernel under a monomial 
embedding. The pulled-back Fubini-Study kernel is very simple to analyze and is the raison d'etre of 
our method. The Toeplitz-Fourier multiplier requires more work. In Proposition |2.2| , we prove that this 
operator is a T m+1 -invariant Toeplitz operator of order m on Xp (modulo a smoothing operator). We apply 



Proposition 2.2 to obtain the goal of this section, Lemma 2.3, which allows us to use the method of stationary 



phase to obtain the mass asymptotics of Theorem g. 

2.1. The exact formula. Our exact formula for the Szego kernel of a toric variety involves two ingredients: 
the first is the kernel 

u^(x, y ) := -LnfVCzWO/)) = (t>p(*),TF(y)) = £ ftfoMto , (47) 

a£P 

where S = #P — 1 and ip : Xp — > S 2S+1 is the lift of the monomial embedding Mp «— ► CP 5 , as described in 



1.3.We shall choose the constants c a = 1 for all a S P, so that = Xa- (Except for Proposition 2.1, all 



our results hold without change for arbitrary {c Q }.) For simplicity of notation, we shall write 

nf(x,») = [nf'(x, y)] N . 

We note that nf (x,x) = 1. 

Second, we introduce two sequences of functions Vn and Qn on NP for each iV = 1, 2, • • • : 
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(i) The 'partition function' V N {a) = #{(/?i, . . . , 0n) ■ 0j € P, [3i H h /3 N = a}, where a £ NP. It 

counts the number of ways of writing a lattice point in NP as a sum of N lattice points each in P. 

(ii) The norming function: Qn{o) := J Xp \x P (x)\ 2 dVolx P (x). 

In our application, we only consider the product 

V N Q N {a) = f f n^(z,e^z)e- i ^>d^dV(z). 

J M P JT 

We observe that each sequence of functions can be re-defined as a single function on the homogenized 
lattice cone Ap = Uv=i NP. By definition, NP C Z" i+1 is obtained by homogenizing a £ NP to a Np := 
(Np — \a\, ai, . . . , ct m ). We then define P, Q : Ap — > R + by: 

7?^) 1=7^ (a), Q(a JV ) = Qjv(a) ■ 

Since the monomials Qjv(a) _ ^Xa form a complete orthonormal basis of H 2 (Xp), they define the eigen- 
values of a special kind of operator on H 2 (Xp). 

Definition: An operator HFU on Ti 2 (Xp) will be called a Toeplitz- Fourier multiplier if it satisfies the 
following (equivalent) conditions: 

• F may be expressed as a function F(D) of the commuting system of operators D = (Si, . . . , S TO , J^). 

• Its eigenfunctions are the monomials x P . 

Thus, we put: 

V{D)x P a =V{a N )x P a , Q(D)x P = Q(a N )x P , a £ NP. (48) 

The following explicit factorization formula exhibits the Szego kernels II^ P as the composition of two 
simple operators. It is the basis for our analysis of the analytic continuation Il^ p . 

Proposition 2.1. We have: 

n™ p = (PQ)- 1 nf . 

Proof. First, we have by definition, 

On the other hand, by definition of the partition function, we also have 

n?{x,y)= J2 V N (a) X P (x)W¥) ■ 

a£NP 

We note that the ./V th power of II* fp gives the sum over the correct set of exponents but does not have the 
correct normalizing coefficients. We need to divide each term by Pjv(a)2Af(a) to adjust the coefficients. 
That is just what the Proposition claims. □ 

Since 

(VQ)- 1 ^ = n^(^Q)- 1 n^nf, 

we see that H^ p is the composition of the Toeplitz-Fourier multiplier F±^ p (PQ)~ 1 Tl^ p with the simple 
kernel 11^. 
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2.2. Toeplitz operator approach to n^f p . Our next step is to prove that the Toeplitz-Fourier multiplier 
njy P CPQ) _1 n^ p is a more familiar kind of operator, namely a Toeplitz operator. 

We begin with some background on symbols. We consider symbols a of the form a(z, Dg), where tx(z, N) 
is a semiclassical symbol. Here, we say a £ S k (M x N) is a semiclassical symbol of order k if 

2-1 

<T N (z)=<T(z,N)=N k ^ aj (z)N-i +r l N (z), with \\D n z r l N {z)\\ < C nl N k ~ l (I > 1), 

3=0 

where aj £ C°°{M). It is a smoothing symbol if ||-D™cr^r(z)|| < C n N~ l for all n, I. The Toeplitz operator 
associated to a symbol a is the operator Ha(z, Dg)U, where Dg denotes the symbol of (the S 1 generator). 
Its symbol is the polyhomogeneous function on the symplectic cone £ = {{x,ra x ) : r > 0} C T*X given by 

oo 
3=0 

Since a commutes with the 5* 1 action, we have LlaTI = ~^2 N TLnctnTIn . 

To simplify notation, we shall write Tlx = n^ p . The goal of this section is to prove that Hn(VQ)~ 1 IIn 
is a semi-classical Toeplitz operator in the following sense: 

Proposition 2.2. There exists a symbol ax of order m with principal symbol equal to 1 and a smoothing 
operator Rx so that 

UnIVQ)- 1 ^ = Il N a N Il N + R N . 



As a corollary of Proposition 2.2, we obtain the following lemma, which is the motivation for our approach 
to II|jvp through nj^ p . It provides an essential ingredient in our derivation of the mass asymptotics. 

Lemma 2.3. There exists a symbol ax £ S m {M 1 N) with principal symbol equal to 1 and a smoothing operator 
Rx such that for t £ T m , we have 

IIjv(x, t ■ x) dVolxp (x) = / aN(x)Hi (x,t ■ x) dVo\x P {x) 

X p J Xp 



Rx{y, x)Il^(x, t ■ y) dVolxp (x) dVo\x P (y). 

x P Jx 



Proof (assuming Proposition 2.i ): Applying 11^ = HxiVQ^N to the identity of Proposition 2.2, we obtain 

n N = {n N a N n N + R N )n? . 

Now let T t denote the translation operator f(x) i— > /(i -1 • x) on C 2 (Xp). Since [Tt, II/v] = 0, we then have 

T t n N = T t n N (a N + R N )n? = n N T t (a N + R N )n? . 

Therefore, 

/ LLyO -1 ■x,x)dVol Xp (x) = TraceT t ITv = Trace II N T t (a N + R N )Ilf 

JXp 

= Trace T t a N U^ + Trace T t i? w nf 

crjv(t _1 -x)nf (r 1 ■ X,x)dVo\ Xp (x) 

Xp 

Rxir 1 ■ x, t/)nf (y, x) dVolxp (x) dVo\ X p (») . 



' X p J Xp 

The conclusion follows by the change of variables x i— > t ■ x. □ 
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2.2.1. Boutet de Monvel-Sjostrand parametrix. In proving Proposition 2.2 and other statements about Toeplitz 
operators, we need some a priori facts about H Mp which follow from the Boutet de Monvel-Sjostrand 
parametrix construction. In this section, we provide the relevant background. 

The Boutet de Monvel-Sjostrand parametrix is a complex oscillatory integral formula for the Szego kernels 
of a strictly convex CR manifold modulo smoothing kernels [BS]. It can be used to obtain asymptotics of 
projection kernels onto holomorphic sections of powers of a positive line bundle over a Kahler manifold 
(see [Zc, BSZ]). We recall here the basic properties of the Szego kernel for (Mp,Lp). Unfortunately, the 



remainder term is not sufficiently precise for the tunneling theory we have in mind. In the next section, we 
give a special construction of the Szego kernels of toric varieties which is sufficiently precise. The construction 
uses the general properties established in this section. 

Now consider any polarized algebraic manifold (M, L), and assume L is very ample. Choose a basis 
{So, . . . , Sd} with d = dim H°(M, L) and let $ denote the associated embedding into projective space. 
That is, we write Sj = jj&h relative to a local frame and put $(z) = [fo(z), . . . , fd(z)]- Recalling that 
L = $*0(1), we can equip L with the metric <E>*/i FS . We also give M the metric lu = $*ui FS , and we let 
denote the orthogonal projection onto H°(M, L N ) with respect to these metrics. Also, let n = n^r 
denote the Szego kernel. 



It was proved by Boutet de Monvel and Sjostrand [BS| (see also the Appendix to [BGJ) that II is a complex 
Fourier integral operator of positive type, 

lie I°{X x X,C) (49) 

associated to a positive canonical relation C. For definitions and notation concerning complex FIO's we refer 
to [MS, BS, |BCfl . The real points of C form the diagonal A^xs in the square of the symplectic cone 

£ = {ra x : r > 0,x E X}, (50) 



where a is the connection form. We refer to [BG] (see Lemma 4.5 of the Appendix). Moreover, in 
parametrix is constructed for II, from which it follows (see Zc , BSZ|) that 



/>00 />27T 

U N {x,y)^N / e N (- ie+t ^ ri>x ' y »s(rex,y,Nt)d6dt, (51) 
Jo Jo 

where s(x,y,t) ~ Y^kLo t m Sk{%, y) € S m (X x X x R + ) is a classical symbol of order m. Here, '~' means 
modulo a rapidly decaying term (i.e., a term whose C J norms are 0(N~ k ) for all j, k). 

To describe the phase in (|5l|), we let be a nonvanishing holomorphic section of L over an open 
U C M, and consider the analytic extension' a(z,w) of a(z,z) := a(z) — ||ei(z)||^ 2 in U x U such that 
a(w, z) — a(z, w) on U x U. Using coordinates (z, 9) for the point x = e ie a(z)^eL(z) € X, we have 



tjj(x 1 ,x 2 ) 



-1 + e 



i(6 1 -6 2 )_ 



a(zi,z 2 



(52) 



y/a(z 1 ) y /a(z 2 ) 

Again assuming that the metric to on M is the pull-back of lu fs , we claim that the phase equals: 

^(X,y) = -1 + (i(x),i(y)) , (53) 
where i : X — » C 2d+1 is the lift of $ given by 

L (z,6)=e* e (^|/ J (z)| 2 )^(/ (z),...,/ d (z)). 

To see this, let us recall the Szego kernel of the hyperplane section bundle 0(1) — > CP d over projective 
space with the Fubini-Study metric. We take U = {z ^ 0} w C d , and we consider the local frame 
e = z . Using the local coordinates [1, z\, . . . , z m ] i— > (zi, . . . , z m ) S C d , we then have a CP (z) = ||e(2;)|| _2 = 
1 + \ z j\ 2 ■ Hence a C¥ has the real-analytic extension 



a Cri (z, W ) = l + £, 



(54) 



3=1 
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and (|53| ) follows. In the case of a toric variety Mp, ( |53|) becomes 

i>(x, y) =~ 1 + Y1 ™^( x )™%{y) ■ 



2.2.2. Complex Fourier integral operators. We now use the theory of complex Fourier integral operators 



| M£ | (or alternatively, Toeplitz operators [BGQ to prove that the Toeplitz F ouri er multiplier II (VQ) 1 II 



is a Toeplitz operator of order m modulo a smoothing operator (Proposition UM. Recall that a Toeplitz 



operator in the sense of [BG] is an operator of the form LL4LI where A £ >L m (X) for some m. Here and 
in the following, ^ m (X) denotes the space of pseudodifferential operators of order m on X . A smoothing 
Toeplitz operator is a smoothing operator of the form Hffll. 

It follows then that IIjv = TatII^ where TV is a Toeplitz operator. Since it is of some independent interest 
and requires no extra work, we prove this as a corollary of a general result valid for any polarized algebraic 
manifold (M,L). We assume L is very ample, and let i : M — > CP d denote the holomorphic embedding 
induced by H°(M,L). Let Ii b {x,y) — U(l(x) , t(y)) denote the pullback to X of the Fubini-Study Szego 
kernel of CP rf . Also, let IL M or more simply II denote the Szego kernel for M with respect to t*w FS . 

About complex Fourier integral operators A (or Toeplitz operators), all we need to know are the following 
basic facts: 

• A £ I" l (X x X,C) possesses a principal symbol a a which is a half-density (times a Maslov factor) 
along the underlying canonical relation C. (In the Toeplitz case, it is a symplectic spinor.) 

• We can compose operators in I™(X x X,C) on the left and right by elements B £ ty k (X) and uab = 
o~ao~b = o~ba- (The same is true of Toeplitz operators.) 

• If A £ I™(X x X,C) and if a A = 0, then A £ I'^ 1 {X x X,C) (and also for Toeplitz operators). 



II and II t are elliptic in that their symbols are nowhere vanishing. See [BS] for the symbol in the 



complex FIO sense and [BG, §11] for the symbol in the Toeplitz sense. 



Lemma 2.4. Let II denote the Szego projector associated to the metric i*w ws . Then, there exist A £ \E' m (A) 
with [A, Dg] = such that LT ~ LL4IIIL = IL4II t modulo smoothing operators ILRII. 

Proof. As mentioned above (fl9|), LT is a complex Fourier integral operator associated to a positive canonical 
relation C, whose real points form the isotropic relation As x s C T*X x T*X (the diagonal). We observe 
that also II t £ I* (X x X, C). This follows immediately from the fact, show by (p3|), that LI and IL are complex 
Fourier integral distributions with precisely the same phase functions. Since the underlying canonical relation 
is parametrized by the phase, they both belong to the same class of Fourier integral distributions. 

Now, the principal symbol an of LT, viewed as a complex Fourier integral distribution, is a nowhere 



vanishing 1/2-density on C which is computed in [BS, Prop. 4.8]. Alternatively, viewed as a Toeplitz operator 
in the sense of [BG| (see Chapter 11), its symbol is an idempotent symplectic spinor . Similarly, the principal 
symbol <rn t of IT, is the pull back under l of the nowhere vanishing symbol of LT CP . 

By our normalization, II t has order — m (since its amplitude is a constant independent of N). We therefore 
begin by seeking Aq £ 1 $ m (X) such that [Aq, Dg] — and such that II — IL4oIL is of order —1. We first find 
oq £ C°°(M) such that on = aooir an d choose Ao so that [Aq, Dg] — and so that <ta = do- Existence of 
such an Aq follows by ellipticity of IL and by averaging; see also ]BG| , Prop. 2.13]. Thus, the principal symbol 
of order of LT— LTAoLT,, equals zero, i.e. II— nAoIL £ I^ 1 (X x X, C). We denote its principal symbol by U-y. 
We then seek A_ x G * m (^) so that [A_i,D e ] = 0, and so that II — IL4 IL - LM-iZ)^ 1 ]!, £ I~ 2 (X x X,C). 
Here, we note that UDgll is an elliptic Toeplitz operator; since nn,, = n,,, the expressions Z?^" 1 n t are well- 
defined. It suffices to choose a_i = o~a_ 1 £ C°°(M) so that a^iari, = c— i- We continue in this way to obtain 
a_j £ C°°{M) always using that o~n L is nowhere vanishing. By a Borel summation argument, we can find 
A £ ^ m (X) with the above commutation properties so that UAH ~ X)j=o ^A-jD^I! is a Toeplitz operator 
of order —M — 1. Then 

n - IL4IL £ r°°{x x X,C). 

□ 
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We remark that an alternative to the observation that the two complex FIO's have the same phase function 
is that, by the choice of the Kahler form i* cu FS , i is a symplectic (as well as holomorphic) embedding M C CP d . 
Hence, the pull back operator l* carries the class of FIO's in the class of n cpd to those in the class of IT. 



2.2.3. Semi-classical Toeplitz operators. To complete the proof of Proposition 2/2, we need to rela te the 
complex FIO IL4IF to semi-classical Toeplitz operators. The relevant analysis already exists in [GuS]. 

We will need the following asymptotic formula for symbols. In the language of Berezin- Toeplitz operators, 
we are essentially computing the Berezin transform between covariant and contravariant symbols of a Toeplitz 
operator (see | RT for background). 

Lemma 2.5. Let cr/v ~ N k Y^jLo S-jN -1 * be a semiclassical symbol of order k. Then there exists a complete 
asymptotic expansion 

i-i 

N- m (U N (T N U N )(z,z) =J2b-A z ) Nk ~ J +r l N (z) (* > 1) . 

3=0 

where bo — so,b-\ — Asq + s_i, . . . , and in general where b^j is a sum of differential operators applied to 
80,8-!,..., s-j. Also, \\Dy N {z)\\ <C nl N k ~ l . 



Proof. Apply the method of stationary phase using (51)-(B2) exactly as in the proof of |Zcj, Theorem 1] 
(which is the case a = 1). □ 

We then have: 

Lemma 2.6. Let (M,L) be a polarized algebraic manifold as above. Then there exists a symbol ajv € 
S~ m {M, N) with principal symbol equal to 1 so that 

nf = n N a N n N + R N where ||-R;v||hs = 0(iV" fe ) Vfc . 

Here, \\Rn\\hs = Trace R* n Rn- 



Proof . By Lemma |2J, there exists A e *~ m (X) such that IT = II All. Since [A,D e ] = 0, it follows by 
HGu2| that there exists a symbol with \\H^ — HncinIInWiis = 0(N~°°). We may determine on by using 
Lemma 2.5. Indeed, we have 

nf (z,z) = 1 - U N a N U N (z,z) (55) 

modulo functions rjy which tend to zero rapidly in C k (M). It follows that o,q = 1 (and then the rest of the 
coefficients may be determined recursively, e.g. a_i = 0, a_2 = —D2C10, and so on). □ 



2.2.4. Completion of the proof of Proposition 2.L In the toric case, IIa^PQ) - IIat is the inverse of 
(V Q)Hn = nf on H 2 (Xp). From Lemma 2.6 we may write IlAr('P<2)Il7v ~ IIn<1n^-n modulo smoothing 
operators. 

We note that we may invert nail in the class of Toeplitz operators; i.e., there exists a symbol <tn such 
that 

TLno-nTIn o HNajyllN ~ n^r (56) 

modulo smoothing operators. Such an inverse symbol exists since 00 = 1. The algebraic formalism in which 
the inverse is calculated is that of *-products of semiclassical symbols. We recall that composition of Toeplitz 
operators defines a *-product on semiclassical symbols by the formula 

UnclnTIn o TlNbiyTlN ~ TIno-n * frjvlljv- (57) 

The formula for * may be worked out directly from the parametrix J5l| ) and the inverse can be computed 
from this formula (see [Gu2, ptfl). 

By (56) we obtain a symbol o~n (with principal symbol equal to 1) such that 

ll N a N U N Il N (VQ)U N ~ U N . 

Multiplying both sides by Un (V S^n^ , we conclude the proof. □ 

Remark: We emphasize that the distinguishing features of the toric case in Proposition 2.2 and Lemma 2.6 
are the exact factorization and the fact that the operator n(7 : 'Q)n mediating between n and n t is invertiblc. 
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This is due to the fact that sections of L generate the ring ®^ =1 H°(M, L ) in the toric case. In general they 
do not and the exact representation II = BTl L in the toric case is only valid modulo smoothing operators. 



2.2.5. Alternate proof of Lemma 2. 6 . One could avoid using the calculus of complex Fourier integral operators 



or Toeplitz operators in the proof of Proposition 2.2 by further developing the calculus of semi-classical 



Toeplitz operators directly from the Boutet de Monvel- Sjostrand parametrix as follows: 

We first simplify the expression (|5l]) by using the complex method of stationary phase to eliminate the 
integrals in the parametrix. The critical point set of the phase $(#, A; x, y) := — iO + Xip{rgx, y) is given by 

d x ${9,X;x,y) = ijj{r e x,y) =0 ^=> e l ° {l{x),l{ v )) = I, 

(58) 

dg<£(9, X;x,y) = -i + d e Xip{r e x,y) =0 <^=> Xe zS (l(x), i(y)) = 1. 

It is easy to see by the Schwartz inequality that a real critical point exists if and only if x = y, in which case 
9 = 0, A = 1 and we obtain the familiar expansion along the diagonal. When i/i/we deform the contour 
to \(\ = e T (C = e l8+T ) so that e te+T (i(x), i(y)) = 1. This is possible as long as (t(x),i(y)) ^ 0, as happens 
near the diagonal, where the parametrix is valid. Because the phase is linear in A it is clear that the critical 
point is non-degenerate if and only if (t(x), t(y)) ^ and that the Hessian determinant equals \(l(x), i(y))\ 2 . 
On the critical set the phase equals —i9 + r = log(6(cc), i(y)), hence we have 

U N (x, y) = e Nlo ^ x ^y»S N (x, y) + W N (x, y), (59) 

where Sn(x, y) ~ SfcLo N m ~ Sk(x, y), and where Wn(x, y) is a smooth uniformly rapidly decaying function. 
Here, we have absorbed the remainder in the parametrix construction as well as the remainder in the 
stationary phase expansion of the parametrix in Wn- Note that the first term may be smaller than the 
second outside a tubular neighborhood of radius N~ x / 2 of the diagonal. 



As above, we use Lemma 2.5 to find a symbol qn so that (55) holds and hence the kernels II ^ and IljvaA/TI 



agree on the diagonal modulo smoothing symbols. Next, we note that both H^(x,y) and Hno.n^n(x, y) 
are complex oscillatory functions with common phase 

y(z,w)=log(i(x),L(y)). (60) 
In the case of , this follows from (p7[). Indeed, we simply have: 

nf (x,y) = e ^i°sW*).%)>. (61) 
In the case of IIjvaArlljv, we apply the method of complex stationary phase to the integral formula 

U N a N Il N {x,y) ~ / e N ^ (u ' x ^S N (x,u)a N (u)S N {u,y)dV(u), (62) 

JMp 

coming from (|59"|), with 

*(u; x, y) = (l(x),l(u)) + (l(u), i{y)). 

It follows that there exists an amplitude An defined near the diagonal such that 

U N a N n N (x,y) = A N (x,y)e Nl °^^ + V N (x,y), (63) 

where Vat is a new smoothing operator. 

Recalling that X c L* , we extend A^ to L* x L* so that it is of the form 

A N (z,X;w,X') = (XX') N A N (z,w) , 

where we use a local holomorphic frame to write x — (z,A), y — (w,X) £ L* . We note that A^{z,w) 
is holomorphic in z and anti-holomorphic in w near the diagonal. To see this, we first conclude from the 



construction in [BS] using the ^-product that the symbol s(x,y,t) ~ /CfeLo *™ ks k( x ,y) m ( pl| ) extends 



to a symbol on L* x L* that is holomorphic in x and anti-holomorphic in y. It follows by the stationary 



phase method described above that the same is true for the symbol SN(x,y) in (59), and then by (|63j) 
that the same is also true for An as claimed. (Note that in terms of coordinates (z, 9) on X, the function 
z i— > An(z, 9; w, 9') — e lN ( e - e ) a(z)a(w) An(z, w) is not holomorphic in z.) 

The amplitude of 11^ — 11^0^11^ becomes (AA') JV (1 — A N (z,w)), where we can write 1 — A N ~ B a + 
BiN^ 1 + B 2 N~ 2 H . Since Ilf - 11^0^11^(2, z) = 0(7V~°°) by our choice of a N above, it follows that 
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Bj(z, z) = (for all j). Since Bj(z, w) vanishes on the diagonal and is holomorphic in z and anti-holomorphic 
in w, it must be identically 0. Hence, — Hjyaj\[H]\[(x,y) is a smoothing operator. □ 

3. Mass asymptotics 

In this section, we prove a precise asymptotic formula for the conditional Szego kernel on the diagonal 
(Proposition 3A ) , which yields the mass asymptotics of Theorem [| 



First, we discuss formula (g) for the expected mass density. Let P denote a Delzant polytope in K m . 
Recalling the definition (||) of the conditional probability measure 7 p |p on the space ff°(CP™ 1 , 0(p), P) of 
polynomials with Newton polytope P, we see that the expected value of the mass density with respect to 
l P \p is given by: 



Since the \ a are independent complex random variables with variance 1 (i.e., E|p(A Q ,A^) = 6@), we have: 

E!p(|/(*) = E%!#= n l^^)- (64) 

aGP IIXa|1 

It then follows by expressing the Gaussian in spherical coordinates that 

v»A\mO = ^e |p (\m\i s ) = ^n |p( z,z) . 

Replacing P with NP, we obtain formula (pf). 



The number of lattice points in the polytope NP is given by the Riemann-Roch formula (see [Fu|): 
#(NP) = dim H° (CP m , O (Np) , NP) = dimH°(M P ,L%) 

m M k 
= x(Mp,L$) = ^deg[ Cl (Lp) fc UTodd m _ fc (Mp)] — 

= Vol{P)N m + • • • + deg Todd m (Afp) . 
Here, we used the Kodaira vanishing theorem and the fact that 

Vol(P) = Vol(Afp) = / —lop = — deg{ Cl (L P ) m } . 
J Mp to! to! 

(This immediately yields Kouchnirenko's Theorem, since deg|ci(L) m ] equals the number of points in the 



intersection of to generic divisors of a very ample line bundle L; see [A.t|.) The holomorphic Euler charac- 
teristic x(Mp, Lp) is also known as the Hilbert polynomial of the polarized manifold (Afp, Lp). Combining 
(^J) and the Riemann-Roch formula, we obtain 

Vu NP (\f(z)\l s ) = ,J tN J ^\np(z,z). (65) 

The mass asymptotics of Theorem || is an immediate consequence of ([35|) and the following asymptotic 
expansion of the conditional Szego kernel on the diagonal: 

Proposition 3.1. Suppose that P is a Delzant polytope in R m . Then: 

i) For z in the classically allowed region Ap, we have 

m 

n lNP (z,z) = l[(Np + j) + R N (z) , \\R N \\cHK)=0(e- x * N ) Vfc, 

for all compact K C Ap, where Xk > 0. 

ii) On each open forbidden region 1Z° F , 

U lNP (z, z) = N ^e- N »W [c F (z) + c({z)N- 1 +■■■ + c F k (z)N- k + R F k {z)} , 
where r = dim/* 1 and 
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(a) 

(b) 

(c) 
(d) 

(e) 



E C°°{n° F ) and c% > on TZ F ; 



l-Rfe lies' (if) 



1 ), for all compact K C 1Z° F and for all j, k; 



= 0(N- k ' 
b > on C* m \ Ap; 
b is given by formula ( \1L V , 
b G C^(C* m ) (with b — on Ap ), and b is C°° on each closed region IZj 



We remark that the asymptotics of Theorem 3.1 are uniform away from the transition points only. To 



take care of the transition points, we shall also prove the following local uniform convergence result on all of 



Proposition 3.2. Let P be a Delzant polytope. Then 



— logH\ NP (z,z) -» -b(z) 



uniformly on all compact subsets of C* r 



This section is devoted to the proofs of Propositions 3.1 and 3.2. We begin with part (ii) of the former, 
which is the major part of the section. The proof begins by using Lemma 2.3 to obtain an oscillatory 
integral formula for the conditional Szego kernel H\np(z, z). We then observe that the integral is rapidly 
decaying when z is in the classically forbidden region, while for z in the allowed region, the phase is of 
positive type with a nondegenerate critical manifold, and hence H\np(z, z) has an asymptotic expansion. 
In §3.2, we consider z in the forbidden region and we seek deformations of the contour of integration so 
that the (analytic continuation of) the phase picks up critical points. In §3.2.1, we formulate necessary 



geometric conditions for the existence of critical points where the phase has maximal real part and we show 
that there is a contour for which these conditions are satisfied. In §§3.2.2- ^.2.5 , we show that the conditions 
are sufficient and, with the added assumption that z is not a transition point, the resulting critical points 
are nondegenerate (Lemma 3.S). We complete the proof of Proposition |3.1|(ii) i n §§3.2.4-3.2.5. In §3.2.5, we 
use 3.1(h) to prove Proposition 3.2, which we then use to prove 3.1 (i) in §3.2.6. 



We start with a simple integral formula for the conditional Szego kernel: 



Lemma 3.3. Fory 1 ,y 2 £ C" 



have 



n\Np{yuV2)= [ f n^(t-y 1 ,y 2 )n^(x,t-x)aNo\ T ™(t)dVo\ Xp (x). 



IXp JT" 

Proof. We consider the lattice sums 

XNp(t) := > * e T" 



aeNP 



We easily see that 



^\np(Vi,V2) = J H"p{t ■ yi,y 2 )XNp{i) dVol T m(t) 



Recalling (|32|), we also have 



x f 



Il™ p (x,t-x)dVolx P (x) = 



t a 

i r P|i 2 \Xa(x)\ 2 dVol Xp (x) = X Np(t) 

Xp aeNP " X P 



E 



The conclusion follows immediately from (|67|)-(|6S|). 



(66) 
(67) 

(68) 
□ 



By Lemmas 2.3 and 3.3, we have for z E C* m , 

1 



H\np{z,z) = 



M P JT 



Wjfp (e** • z, z)n* fp (w, • w) d<p d\o\ Mp H 



(27T) 

K N (z) + S N (z) , 



(69) 
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where 



K N (z) = 
S N (z) = 



(2tt)' 
1 



M P JT 



If^p (e** ■ z, z)a N (w)U^ N (w, e l * ■ w) dip dw 



(70) 



Mp J M P JT 



n^ r p " (e^ • z, z)R N (w, ?y)rif N {rj, e** ■ w) dtp dw d v , (71) 



where <7jv, Rn are as in Lemma |2.3| , and n^ /jv is given by (f47j). (The above are really functions of x — (z,9), 
but since they depend only on z, we dropped the variable 9. In (71), we replaced integrals over Xp with 
integrals over M P by setting R N {w,rj) = jf* R N ((w, ft), (77, 2 ))e lJV ( ei - e2 ^ff .) 
By (40), we have 



J Q I<P 



Recalling (}47|), we can then rewrite (|70|): 

y ' (27r) m 



M P JT 



(Np)\ 

e N ^' w '^a N (w)dipdVol M p{w) , 



where the phase is given by 

V(<p,w;z) = log ^ e-^ a ^\m?(w)\ 2 + log ^ e*< a ^ |toP s (z)| ; 



(72) 
(73) 



aeP 



\a\<p 



and cin(w) = ^^p"^ &n{w) is a symbol of order 2m. 

We fix z £ C* m . It follows from the triangle inequality and (Bq) that 



^(p,w;z) 



log|E Qe pe-^>KH| 2 |+log Ela 



|a|<p' 



(74) 



< 3?*(0,w;z) = , 

with equality if and only if (a, </?) = (0, = for all a G pS, or equivalently, if ip = 0. This says that 'J is a 
phase function of positive type, and w) = if and only if <y9 = 0. Hence we need consider only critical 

points with ip = 0. 

Suppose that u> = e* e e p S C* m . We note that \& is independent of 9; hence the critical set is T m -invariant. 
So we must find the set of w where 



We have 



D v V(Q,w) = Q, D p *(0,w)=0. 
E aeP e- l{a '* ] \m p a {w)\ 2 a E agpS ^ ] K s (z)| 2 a 



Hence the D v equation gives 



- (ip(w) = fis(z) ■ 
p 



(75) 



(76) 



Note that 

D p V(0,w) = D p = 

automatically by (pOh. However, we will need to compute D p ^. To do this, we note that 



aGP 



E 



a£P ' 



m a (w)p 



Recalling that to = e l6 e p , we then obtain 



y 



eP' 



,-»(a,p) \-P 



fh^{w)\ a 



/xp(ui) 



(77) 
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Let C% C T m x Mp denote the set of critical points of ^ where ?RA> = 0. By @ and @, we have: 

Cl = {<p = 0,we[i P 1 ( t iz(z))}^T m . (78) 
Equation ( |78| ) leads to the important conclusion: 

There are no critical points with real part equal to unless z S (/is)~ 1 (^-P)- 

As we stated in the introduction, we refer to the region Ap = {pz,)~ 1 (^P°) as the classically allowed region, 
and its complement in C* m as the classically forbidden region. We now break up our problem into these two 
regions. 

3.1. The classically allowed region. In this section, we illustrate our approach by showing that 

n lNP (z,z)~N m [c (z) + c 1 (z)N- 1 +c 2 (z)N- 2 + ■■■], for z E A P . (79) 

In §3.2, we will derive our exponentially decaying, asymptotic expansion of H\ NP (z, z) on the forbidden 
region, which we will use to obtain the precise formula of Proposition 3.1 (i) on the allowed region. (Although 
we do not need ( |79| ) in the derivation of our precise formula in §3.2, we shall use the methods and Lemma 3.4 
from this section in our proof.) 

Suppose that z e Ap = {ps)^ 1 (^P°) is in the classically allowed region. Since the phase function is of 
positive type and equals zero on the critical set, in order to apply the method of stationary phase it suffices 
to show that the Hessian of ^(ip,p) is nondegenerate on the critical set. We write the Hessian as a block 
matrix: 

A B 



where 

d 2 ^ d 2 ^ <9 2 * 

^ = ^nr(°^)' B jk = ——(o, P ), c jk = ——(o, P ). 

oipjdipk oipjdpk opjdpk 

Differentiating ( |75| ) and recalling that pp(w) — J^aeP \™a( w )\ 2a i we obtain 

A = fip(w) <g) pp(w) - ^2 \™a( w )\ 2a ® a + Mps( z ) ® P P T,{z) - ^2 I^S S ( z )| 2q; ® a ■ 

a£P \a\<p 

(For a vector v — (vi, . . . , v m ) € R m , we identify v (g> v with the symmetric matrix (VjVk)i<j,k<m-) 
Thus on C%, we have by @, 

A = 2p 2 pt,{z) ® ^{z) - \ fh^{w)\ 2 a®a- ^ \rh p ^ {z)\ 2 a ® a . (80) 

a£P \a\<p 

On the other hand, differentiating (0), we obtain 

(81) 



B = —2i ^2 \™a( w )\ 2a ®ot — Hp(w) ® pp(w) 

\aeP / 



Finally, since ^(0,p) = 0, we have 

C = . 

Therefore 

detL> 2 *| (0!p) = (detB) 2 . (82) 

In order to apply the method of stationary phase to conclude that the integral (72) admits a complete 
asymptotic expansion in powers of N^ 1 , we must show that det D 2 ^\t p \ ^ 0. This follows from ( [32"| ) and 
the following fact: 

Lemma 3.4. The real symmetric matrix 

S w := 22 \™ot ( w )\ 2a ® a ~ Mp( w ) ® Hp(w) 

aeP 

is strictly positive definite for all w € C* m . Furthermore, if fip(w) lies in a face F of P, then S w is 
semi-positive and its eigenspace is Tp~ ; in particular, S w is positive definite on Tp. 
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Proof. We must show that (A <E> A, S w ) > for A E (M m )' \ {0}. Consider the vectors u,v € M #p given by 
u a = \m£(w)\, v a = \m%(w)\\(a). Since ||u|| 2 = J2aeP \™-u( w )\ 2 = 1 an d ^p(w) = J2aeP ( w )I 2q; ' we 
have 

(\®\,S W ) = I^M| 2 A(a) 2 - (j2 \™« H| 2 M«)) 

= |H| 2 -(u,v) 2 = \\u\\ 2 \\v\\ 2 - (u,v) 2 > 0. (83) 

If w € C* m , then the Cauchy-Schwartz inequality in ( |S3|) is strict since |m p (?z;)| ^ for all a <G P D Z 
and A(a) is not constant on P n Z. On the other hand, if /J,p(w) lies in a face F, then by Proposition 
l^aC 10 )! / O'SaEf, and hence 

(A ® A, S M ) = <^ X\ F = const. <^ A _L T F . 

□ 

The asymptotic expansion (179) for z in the classically allowed region Ap now follows from the method 



of stationary phase (| Ho, Theorem 7.7.5]) applied to the complex oscillatory integral fl72) and the fact that 
Sn{z) is rapidly decaying. To be precise, the asymptotics of (|7^ ) are determined by the component of the 
critical point set of the phase on which Sft'I' is maximal. Note that 3?\1/ and are constant on components 
of the critical set, and that on its maximal component our phase satisfies 3?^ = 0. The critical point set 
is a manifold and we have just shown that the phase is non-degenerate in the normal directions. Although 



([Ho, Theorem 7.7.5]) assumes the critical point with 3?^ = is isolated, the proof in [Ho] can be modified, 
precisely as in the real case, to apply to phases with non-degenerate critical manifolds. 

Hence the expansion follows from the complex stationary phase method for phases of positive type with 
non-degenerate critical manifolds. Since the critical manifold is of (real) codimension 2m and the amplitude 
djv in (j7^) is or order 2m, it follows that the leading term of the expansion contains N 2m ~ m — N m , as 
claimed. 

3.2. The classically forbidden region. For z in this region of C* m , there are no critical points of the 
phase with real part 0, and the integral is rapidly decaying by (f74|). Perhaps surprisingly, this implies that 
dd\ogTl\Np(z, z) will contribute here to the distribution of zeros. 

To determine the exponential decay rate of the integral (f72|), we shall deform the contour of integration 
to pick up critical points with maximal real part along the contour. To accomplish this, we first complexify 
the real torus T m to C* m with variables Q = e Tj+tlfj . The integral ( |72] ) may be written in terms of these 
variables as 

Kn(z) = j^— I I e N ^^a N (w)f[^idVol Mp (w), (84) 



where 



y c (cw;z) = io g £c Q KMi 2 +iog E n^fwi 2 

a£P \a\<p 

= log J2 e-^ T+ ^\m^(w)\ 2 +log Y, e <a ' T+ ^|mf (z)| 2 . (85) 

aeP \a\<p 



Since the integrand is holomorphic in £ 6 C* m , we can deform the contours in ( J84| ) and instead integrate 
over 

/ / ■ (86) 

JMp ^|Cl|=Tl,- ,\Cm\=T m 

We therefore look for complex critical points of ^c- We note that as before, 

m c (e T+lv ,w) < m c (e\w) , for <p ^ . (87) 
Thus the critical points W with maximal real values on the contour (pq) are those with ip = 0. 
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We now write <&c(t, <p, w) — ^c(e T+lip , w; z). The complexification of J75| ) becomes 

Eg g P e-<°-^> |m^H| 2 a E agpS e<^+^> |m^(z)| 2 a 
S a£pe -(«^)|m^WP E QepS e<<^+^>|m£ S (z)| 2 
E« £ p e-^**") |m Q (u,)| 2 a £ Q£pS e< Q,T+iv) |m a (^)| 2 a 



E Qe pe- <a ' T+lv> |m Q ( W )| 2 Eaepse^+^ha^)! 2 
E ag pe"' <a ^|mg(e-/ 2 - W )| 2 a E agpS ^) |m^( e -/ 2 ■ z)| 2 a 
E Qe p^ ,;<Q ^l^(e-/ 2 - W )| 2 SQepS e»<a, v >|aS E ( e -/2. z )|2 ' 

for to € C* m . Hence by continuity, we see that 

. D = E ag P^ <a - y> l^(e- /2 -^)| 2 a _ E agpS ^)|m^(e-/ 2 -z)| 2 a 

* V C E Qe pe- l ^|m^(e-/ 2 -u;)| 2 e ^> |^ (e r/2 . z) |2 ' [ > 

for all w € Mp. In particular, 

iD v #c(T, 0, «j) = /i P (e- T / 2 • w) - M P s(e T/2 • z) . (90) 
Hence the equation D^ciT, 0,w) — is equivalent to 

-Me- T,2 -w)=^{e T l 2 -z). (91) 
P 

To compute the w derivative, we first consider the case where w = e p+ld € C* m . Recalling that 
V e-<°.T+*P>m'(vA\* E a£ pe-^|m a (e-/ 2 . W )| 2 _ ^ P e^\m a (e^/^)\ 2 

we obtain 



_ £q e pe- i<a - y> l™£(e- T/2 -w)l 
2 ^ P ^ E Qe pe- l ^l™^(e-^ 2 -w)l 2 
Setting 95 = 0, 

D p ^ c {t, 0, to) = MP (e- T / 2 u;) - fi P {w) . (93) 
Hence the equation D p ^c\<p=o = gives 

fi P (e- T/2 ■ to) = hp{w) . (94) 

We must show that ( j94[) also holds for critical points to ^ C* m . (These are points with fip(w) £ dP.) To 
do this we introduce the function 

H P (r, W ) := ^y 6 ^"^ 1 ^^ 12 = E e- <a ' T> KH| 2 ■ (95) 

Note that Hp(r, •) is a well-defined function on Mp that factors through the moment map /ip. We also have 

~l + \\e- T / 2 ■ z\\ 



1+ U 



(96) 



H 1 < e (t,z) = Hv(t,zY' 
We have 

* c (r,0,w;z) =logff P (r,w)+logff pS (-T,z) , (97) 

and hence 

D w ^ c (t,0,w;z) = D w log H p (t,w) . 

Furthermore the vector fields 

^ : =4 = ^Ks|-) (98) 
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extend to Mp as the generators of the M™ action (induced by the C* m action). Hence (93) is equivalent to 
the equation 

{X u . . . ,X m )logH P {T,w) = Li P (e- T / 2 w) - vpM, (99) 

which is invariantly defined on all of Mp. 

Thus we have seen that critical points of i?p(r, •) must satisfy (Q). We shall show later that (|94|) is also 
a sufficient condition for w to be a critical point of Hp(r, •), and hence (|9l]) and (|94| ) give necessary and 
sufficient conditions for (0,w) to be a critical point of ^c( T , •)■ 

The following fact is useful in describing the critical set. 

Lemma 3.5. Let r G R m , w G Mp. Then the following are equivalent: 

a) /ip(e r • w) = fip(w); 

b) e rr • w — w, for all r G K; 

c) r is perpendicular to Tpm, where F w is the face of P containing fip(w). 



Proof. Let w = (w a ) ae p G Mp C CP# P 1 . By Proposition 1.1, we see that w a = if and only if a ^ F w , 
and thus for any t G C* m , 

t-w = w 4^ 3c such that t a w a =cw a Va£ P & t a = c Va G n Z m . (100) 



Suppose (a) holds. Then there exists e lip G T m such that e T+lv ■ w = w. Hence by ( ji00| ), e^ +i * , ' a > = c and 
thus e< rr < Q > = |c| r , for all a Gf*nZ m . Therefore (b) holds. 
Furthermore, it follows from (100) that (b) is equivalent to 

(r,a-a') = Va,a'eP»nZ ra , (101) 

which is a restatement of (c). □ 

Remark: Let r € M m . It follows from Lemma |3.5| that the set A T of points in Mp satisfying (Q) consists of 
the inverse image under /ip of the union of faces of P that are orthogona l to r. The set A T is a (possibly 
empty) smooth algebraic subvariety of Mp. It then follows from (^) and (101) that Hp(r, •) is constant on 
each connected component of A T . Indeed, on a component A of A T , we have Hp(r, •) = e - ( Q > T \ where a 
is a lattice point in (ip(A°). 



3.2.1. The normal bundle. Although equations (91 ) and ( |94| ) are necessary and sufficient conditions for (0, w) 
to be a critical point of the phase on the r-contour (as we shall show), they do not guarantee that the real 
part of the phase is maximal at that point. To obtain maximality, we need to strengthen condition ( |94| ) to 
require that — r lies in the 'normal bundle' of P at w z . 

Recall that the normal cone Gp of a face F of P is given by 

C F = {ue K m : (u,x) = sup(u,y), Vx G F} . (102) 

y£P 

We now consider compact polytopes in R m with arbitrary vertices that are not necessarily integral or even 
rational. We define the faces of such polytopes as in the integral case, and we use (102) to define the normal 
cones. 

Definition: The normal bundle Af(Q) of a (not necessarily integral) convex polytope Q C K m is the subset 
ofTgrn = K m x R m consisting of pairs (x,v), where x G Q and v is in the normal cone Cf of the face F of 
Q that contains x. (Note that Af(Q) is not a fiber bundle over Q.) 

The normal bundle Af(Q) is a piecewise smooth submanifold of Tr™ ; it is homeomorphic to R m via the 
'exponential map' 

E Q : Af(Q) -► R m , £ Q (x,v)=x + v, x g F, c g C f (F a face of Q) . 

It is easily seen that £q is a homeomorphism and is a C°° (in fact, linear) diffeomorphism on each of the 
'pieces' F xC F C J\f(Q). 

We shall use the following elementary, but not so well known, fact about the invertibility of Lipschitz 
maps. We let B £ (x ) = {x G M m : \x — x \ < e} denote the e-ball about x G R m . 
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Lemma 3.6. [Fan] Let f : U — > R' m be a Lipschitz map, where U is open in R m . Then f has a local 
orientation-preserving Lipschitz inverse with Lipschitz constant L at a point xo € U if and only if there 
exists e > such that 

i) liminf„_>o \fi x +v) — f(x)\/\v\ > 1/L for all x G B £ (xq), 

ii) det f'(x) > for all x G B £ (xq) such that f is differentiate at x, 

iii) f(xo) £ f(dB s (x )) and deg [/ : OB £ (x ) - R m \ {/(a*)}] = 1- 

By the deg ree in (iii), we mean the degree of f*(dB e (xo)) G i? m _i(R m \ {/(xo)},Z). The hypotheses (i) 
and (iii) in [ Fan| differ slightly from those above, but Fan's proof of sufficiency uses only (i)-(iii). (Necessity 



is obvious.) Hypothesis (i) is given as a lemma in [Fan]. Hypothesis (iii) above is replaced in [Fan by the 



condition that the index of / at xo is 1, which is easily seen to be equivalent to (iii) under the assumptions 
(i) and (ii). 
We now let 

7V° =Af(±P) n (E° x R m ) . 



Lemma 3.7. For each point z G C* m , there exists a unique (x,v) G Af° such that 

H E (e v ' 2 ■ z) = x . 



Proof. We note that the R™ action on C* m descends via the moment map to an R™ action on E°, and 
we consider the map 

$ : Af° — > E° , v) = /is {e v/2 ■ x) . 

It suffices to show that $ is a bijection. 

Let £ : E° — > R m be the diffcomorphism given by 

""■(^ "i^;)' (103) 

so that 

Lo^{z) = (log|z!| 2 ,... ,log|z m | 2 ) . 

Thus, writing x = ^(z) we have 

£o$(x,v) = Co nz(e v/2 ■ z) 

= (v 1 +log|zi| 2 ,... , i> m + log|z m | 2 ) 

= v + £(x). (104) 

We first observe that £ o $ is proper: suppose on the contrary that the sequence {(x n ,v n )} is unbounded 
in J\f°, but £ o <&(x n ,v n ) — > a G R m . By passing to a subsequence, we can assume that x n — > a; G ^P. 
Then x° G 9E, since otherwise u" — > a — £(x°). Write u" = r n u n , where r n > 0, \u n \ = 1. We can assume 
without loss of generality that u n — ► u°. We consider the case where Ex° < 1, i.e., /i^ 1 (a; ) is not in the 
hyperplane at infinity. If a; = 0, then u° > (1 < j < to), since otherwise £ o <fr(x n , w") would diverge. But 
u° G C{ }; hence |P C {x : (u°,x) < 0} C R m \E°, a contradiction. Now suppose that x" = ■ ■ ■ = x° k = 0, 
X[ > for k < I < to. Then we conclude as before that u° > for 1 < j < k, and = for k < I < to, and 
we again obtain a contradiction. Finally, if Ea;^ = 1, we can change coordinates (permute the homogeneous 
coordinates in CP m ) to reduce to the previous case. 

Let S = Sip : Af(-P) ^ R m , and let U = S(Af°) C R m . We consider the map / : U -> R m given by 
f o S\j\f° — £ o $, i.e. by the commutative diagram: 

W° E° 
[/ R m 

Since £ o $ is a proper map, / is also proper. Hence to show that $ is a bijection, it suffices to show that / 
is a local homeomorphism and is therefore a (global) homeomorphism. 
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To describe the map /, for each i e E°, we let q x denote the closest point in ~P to x. We note that 
q x S S° ; if x ip, then q x £ d(^P); if x S ^P, then = a;. Furthermore, = (g^, a; — q x ) and hence 

/(x) = x - q x + C(q x ) . 

We shall show that / satisfies the hypotheses of Lemma [D]. Let F be an arbitrary face of ^P n S°. To 
verify (i) and (ii) , it suffices to show that det Df > on the (noncompact) polyhedron 

U F := £{F xC F )f\U = £[{F n £°) x C F ] . 

To compute the determinant, we let T F C W n , N F — Tp~ C M. m denote the tangent and normal spaces, 
respectively, of F. Let x° £ F be fixed. For y £ U F , we have 

y = x + v i — > £(2;) + f , a; — x° € Tp, f G ATp . 

Choose orthonormal bases {li, . . . , Y r }, {Y r+ i, . . . , Y m } of T F , N F , respectively. We let Df denote the 
matrix of the derivative (f\u F )' with respect to the basis {Yi, . . . , Y m } of R m . We have: 

D/=( W 'f!), (105) 



where T is the m x r matrix [Yi • • • F r ]. We have by ( |103| ), 

1 1 m 

(£'(x)l = — + <5*— , i o = 1-Vi j >0, (106) 
V /jfe x J xj ^ 

for x £ E°. Hence £{x) is a positive definite symmetric matrix, it being the sum of a semipositive matrix 
(all of whose entries are and a positive definite diagonal matrix. Therefore, T*£'(x)T is positive definite, 
and hence det Df(x) — det(T* C ' (x)T) > 0, completing the proof that hypotheses (i) and (ii) of Lemma 3.6 
are satisfied. 

Note that (106) implies that the eigenvalues of C'{x) are > 1 for x £ S° and hence by (105), Df{x) is a 
diagonalizable matrix whose eigenvalues are real and > 1. 

We verify (iii) by a homotopy argument: Choose a point q° £ (ip)°. We contract ^P to q°; i.e., for 
< t < 1, we let 

Qt = (1 - t){q } + -P , 
P 

so that Qo = {q }, Qi — ~P- For < t < 1, Qt C S°, and hence we have a map 

$t : JV(Qt) -> S° . 
For < t < 1, we define ft : W n — > K m by the commutative diagram: 

yV(Q t ) E° 

As before we have 

f t (x)=x-ql+L(ql), (107) 

where q l x is the (unique) point of Q t closest to x. The above argument shows that the maps f t also satisfy 
(i) and (ii) of Lemma |3.6| 
We write 

F : (R m x [0, 1)) U (U x {1}) -> E m , (x,i) ^ / 4 (x) , 
where /1 = / : U — > K™. One easily sees that 

ql - q£' <\x- x'\ + |t- t'| , 

and hence F is continuous. Furthermore, F is uniformly continuous on W n x [0,<], for each t < 1. 
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Let H denote the set of t G [0, 1) such that ft ■ K m — ► lR m is an orientation preserving homeomorphism. 
We note that G H since 

fo{x) = x — q + C(q) = x + constant . 

Note that ( |105| ) also holds for the maps ft and hence as above, the eigenvalues of Dft are real and > 1. 
Hence f t satisfies conditions (i) and (ii) of Lemma 3.6 with L = 1, for all t G [0, 1]. 

First we show that H is open in [0, 1). Suppose to G H. Then by Lemma [T^, Lip(/ t ^ 1 ) < 1 and hence 

\x- Xq\ = 1 \f to {x) - ft {x )\ > 1 • 
By the uniform continuity of F on R m x [0, ^^], we can choose e > such that for all x G K m , we have 

\ft(x) - ft{x )\ > \ , for x G dB 1 (x ), \t - t \ < e . (108) 
To simplify notation, we shall write 

deg(/,x ,r) := deg [/ : dB r (x ) - W n \ {/(a*)}] • 



We conclude from (108) that for \t — to\ < e, 

deg(/ t ,x , 1) = deg(/ to ,a;o, 1)) = 1 , 



and hence ft is a local homeomorphism at xo, by Lemma 3.6. By the very first part of the argument 
that / : U — > M. m is proper, we easily see that ft is proper. Since xo G M. m is arbitrary, it follows that 
ft : K m — > R m is a covering map and therefore is a homeomorphism. 

Next we show that H is closed in [0, 1) a nd hence H = [0, 1). Let t n G H such that t n — > to G [0, 1). 
Since /t satisfies condition (i) of Lemma |3j| we can choose e > so that ft {xo) & ft (dB E (xo)). Then for 
n sufficiently large, 

deg(/ to ,a;o,e) = deg(/ t „, x , e) = 1 • 

It follows as above that ft : R m — > M m is a homeomorphism. 

We have shown that ft is a homeomorphism for < t < 1. To complete the proof of the lemma, we must 
show that / = fi is a local homeomorphism. So we let xo G U be arbitrary, and we choose e > such that 
B e (xo) C [/ and /(xo) G" f(dB £ (xo)). Then for t < 1 sufficiently close to 1, we have as before 

deg{f,x ,e) =deg(f t ,x ,e) = 1 . 

Thus by Lemma p.q, / is a local homeomorphism. □ 



As a consequence of Lemmas 3^ and 3/7, we immediately obtain: 

Lemma 3.8. For each z G C* m , there exists t z G E m , w z G Mp so i/iai 

• w(e Tz/2 ■ z) = ~[i P (w z ); 

• (^fip(w z ),—Tg) G J\f°. 

Furthermore, t z , w z are unique ( modulo the T m action on Mp ) and also satisfy §9\ ) and 

We shall show in §§ |3.2.2 |?.2.3| that along the contour given by r = r z , the phase *c has critical points 
and maximal real part at (0,w; 2 ). 

We recall the concept of transition points mentioned in the introduction: 

Definition: A point z G C* m is said to be a transition point if —t z is in the boundary of the normal cone of 
P at jip{w z ), or equivalently if z is in the boundary of one of the regions TZf described in the introduction. 



Remark: It follows from (104) that 2 Log (TZf) = £(F) + Cf, where Log (z) = (log \z\\, . . . , log \z m \). Hence 
we can decompose K m as the disjoint union of the sets C(F) + Cf. If z is a transition point, then 2Log (z) 
must lie in the common boundary of at least two of the sets C(F) + Cf- Figure || below shows the transition 
points in log coordinates (as solid lines) for the case where P is the polytope of Figure ||. 

Note that although Cp' is orthogonal to F' in Figure ||, Cpi is not orthogonal to C(F'). Similarly, in 
Lemma 3.S, if (Ap(w z ) G F' , then the orbit n^(e rTz ■ z) is not orthogonal to F' . 
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Figure 4. Transition points and regions in log coordinates 



3.2.2. Computation of the Hessian. We consider 

C z := {e ie ■ w z : e ie e T m } = {e w ■ w z : 9 e T F } C M P , 

which is a (totally real) submanifold of Mp of dimension r. We are going to show that {0} x C z C T m x Mp 
is the set of critical points of \&c with maximal real part along the t z contour. In order to apply the method 
of stationary phase as in the case where z is in the classically allowable region, one would also like to show 
that {0} x C z is a nondegenerate critical submanifold; i.e., the normal Hessian of the phase is nondegenerate 
along {0} x C z . But the normal Hessian turns out to be degenerate whenever z is a transition point. For 
example, if z € p^ 1 {F) 1 where F is a codimension r face of ^P, then the corresponding normal Hessian 
has an r-dimensional nullspace (as our computations will show) and hence the normal Hessian is degenerate 
whenever /is(^) € d(^P). However, we shall show: 

Lemma 3.9. Let z € C* m . Then {0} x C z is the (unique) component of the critical set of x &c( T z J •, ■) where 
9M/c attains its maximum on T m x Mp. If z is not a transition point, then {0} x C z is a nondegenerate 
critical submanifold. 

Recall that in the previous section, we verified the proposition for the case where z is in the classically 
allowed region. So we now fix a point z in the classically forbidden region. Let F z be the face of P containing 
q z := fip(w z ), and let r denote the dimension of F z . Note that F — q z an open subset of the tangent space 
T F of F C R m . 

Let Y = ^ip (Fz)t which is an r-dimensional complex submanifold, and thus has real dimension 2r. We 
let V C Tm p ,w z be the tangent space to the M™ orbit of w z . Recalling (|98|), we see that 



V 



span{Xi(u; z ), . . .,X m (w z )} = a,jXj(w s ) : (oi, . . . ,a m ) G T F \ 



and hence dim V = r. (Note that V — JT V , where T v is the tangent space to the fiber of fip.) 
The normal bundle N to {0} x C z can be decomposed as follows: 

N = T T m e V®N Y - (109) 

To prove Lemma 3.9, we need the following: 

Lemma 3.10. The normal Hessian of^c{ T z, •) to {0} x C z with respect to a basis for N giving the decom- 
position ( \10!\ ) is of the form 

A 

D^ c (t z ,0,w z 
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where 

i) A is an m x m negative- definite, real symmetric matrix; 

ii) B is an m x r pure imaginary matrix of rank r; 

iii) C is a (2m — 2r) x (2m — 2r) real symmetric matrix, which is negative definite if z is not a transition 
point. 

Proof. Let 

I A B -k 
D^ c (t z ,0,w z ) = B l *" *' 
\ *' *" C 
To simplify notation, we now write r = t z . 

(i) To compute A, we see from ( [B9|) that 

iD^d^a = V" |m^(e~ r/2 • u;)| 2 a + a function of z , (110) 

aeP 

As in our previous computation, we then obtain 

A = /j P (e~ r/2 • w) ® np(e~ T/2 ■ w) - |"i^(e~ r/2 • w)| 2 a<8> a 

By Lemma 3.4, A is the sum of a negative semi-definite term and a strictly negative term, and hence is 
negative definite. 

(ii) To compute B, we similarly obtain 

iD p D^ c \ v =o = Hp{e- T/2 ■ w) <g> Me~ T/2 ■ w) - ^ |ro£ (e- r / 2 • u;)| 2 a ® a , (111) 

where D p = (d/dpi, d/dp m ) = (X 1 , . . . , X m ). By Lemma [TJ, 

(iDp-D^cl^o, A <8> A) > 0, for A G Tp \ {0} , 

and hence B has rank r. 

(iii) To determine C, we need to change variables. Let us recall that Mp is covered by affine coordinate 
charts centered at the vertices of P described as follows: Let v = a be a vertex of P. Since P is simplicial, 
there are exactly m edges (1-dimensional faces) incident to v. Choose a 1 ,... ,a m on the respective edges. 
Let /3 J = a? — a°; then {/3 1 , . . . , f3 m } form a basis for R m . Since we are assuming that P is Delzant, we can 
choose the a? such that {/3 1 , . . . , /3 m } generate T. m . We consider the m x m matrix of integers 

A = (A jk ) = (Pi) . 

Then det A = ±1 and the matrix T := A -1 has integer entries. Consider the (Zariski) neighborhood of v: 
U v := {w G M P : X v(w) ^ 0} = |J{F : F is a face of P, u G P} . 

(The U v cover Mp.) Let 77^ = , J = 1, . . . , m. We note that the monomial mapping 

n = ( m ,...,7 lm ):C* m ^C* m 

is bijective, its inverse given by 

r 3i r, m 

Let F 1 , . . . , F m denote the facets of P incident to v, indexed so that atj £ F- 7 . Let x G F? be arbitrary. 
We then have 

Xai(w) Xai{x) . . 

»7j = -> ^ = as fi P (w) -> x G P J . 

(Since a G -F-?, x Q o(x) 7^ 0.) Therefore 77 extends to an isomorphism 77 : (Uv,v) (C" l ,0). Furthermore, 
the facets incident to -u are the images of the divisors {r/j = 0} under the moment map. (This shows that 
Mp is nonsingular if the Delzant condition is satisfied.) 
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We note that the standard unit vectors e J = r/3 J = Ta^ — Tv. Consider the nonhomogeneous linear map 

f : Z m -► Z m , u i — ► Yu — Tv . 

Then f (a ) = 0, f (a 3 ) = e J (1 < j < to). Consider the polytope Q = f (P). Under the map f, the closed 
facets Fi correspond to polytopes in the coordinate hyperplanes {xj = 0}. Hence the polytope Q lies in the 
upper quadrant {xj > 0, 1 < j < to}. We also let 

A = f - 1 : Z m -> Z m , u i-> Au + v . 

Now let us rewrite the phase in terms of the rj coordinates. We first note that 



Hence, 



H p (t,w) 



E 



<Q < r> |c Q x r(Q ^)l 2 E 



EaeP |CaXr(a)(^)| S 



E 



7€Q I A(7 



) X 7 ( ? 7)I : 



(112) 



We first consider the case where w z — v; i.e., r — 0. Expanding (112) in powers of rjj, we have 



-<u,r) 



H p (t,w) 



-{ex>,r) 



c ,i 2 +e™i 



2 m 



-(•u,t) 



i%i 



3=1 



where c£ — \c a j/c a o\ 2 > 0. Since — r is in the normal cone at i>,we have (—t,v) > (— r, or,-) and hence 
(flj, T ) > 0, for 1 < j < m. Now suppose that z is not a transition point; i.e., — r is not in the 
boundary of the normal cone at {u}. Then r is not perpendicular to any of the rays along the edges incident 
to v; so (f3j,r) ^ 0, and therefore (/3/,r) > 0. Hence 



log H p (t, w)=- (v, t) + ^2 b 3 \Vj I 

3=1 



4(< 



1 < 



This leads to two conclusions: (i) d\ogHp\ w=Wz = and hence (0,w z ) is a critical point of &c(t, ',')', 
(ii) C = Dl log F P | w — ^ 2 is a diagonal matrix with negative eigenvalues h\, b\ 1 . . . , fr^, Note that C is 
singular whenever z is in the boundary of the normal polyhedron of v. 

Now consider the case 1 < r < to — 1. (If r = to, then z is in the allowable region.) We can assume 
without loss of generality that Y — fip{F z ) is given by 



7]i = ■■■ = ?7„ 



= 0, r? m -r+l 7^ 0, 



(113) 



In order to decompose the sum in (112), we consider the sets of lattice points 



E; 



{7 e Q : 7 fe = S{ for 1 < k < m - r} = Q n {e' + Ef= m -,+i ^} , 



for 1 < j < m — r. (These arc the indices of the only terms that contribute to the normal Hessian at w z .) 
We also let E = Qn ({0 m _ r } x Z r ). Since {ft, r) = for I > to - r, it follows that 



and similarly 
We write 



(A( 7 ),r) = (A(e»'),r) = <o»",T) 
(A( 7 ),r) = <A(0), T)=<«,r) 



V = 



Let -E 1 '' c Z m , for j = 0, 1, . . . , m — r, be given by 



{0 m -r} x Eq — E , {0 m _ r } x E" — Ej 



V 7 e Ej , 

V 7 e E a . 

yjm— r+li • • • j ?7m) • 

e-? (1 ^ j '5= ^ — r) • 
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Thus we rewrite (112): 
H p (t,w) = 



^ T) E- ( eE \c Mj) x,(v)\ 2 +£7=T 


e-^' r >E 7eB J^ (7) x 7 (i7)|^ 


H 


E 7S is ic A(7) x 7 wi 2 + Ejrr 


E T eB j |ca( 7 )X7WI 2 


H 



-(f,T) 



E, eBir ^|x,W')l 2 + Er=7 


E^-^lx^")! 2 " 


e" 


^te| 2 + --- 




E^^Mr/oK 


M 2 + - 



wo 

t WO 



(e-^V>-l) |%f + 0(||r/|| 4 ) 



(114) 



where Cj„ € 



and 



WO = E c 'oMv")\ 2 = c' j0 + E c 'jMv")\ 2 > (0 < i < m - r) 



(Note that (p3) implies that H p (t,w) = for w e F z , as we observed earlier, and that (0,w z ) 

is a critical point of *f?c-) Suppose as before that z is not a transition point; i.e., — r is not in the 

boundary of the normal cone at F z . Then ((3j, r) ^ and in fact ((3j, r) > 0, for 1 < j < m — r, since — r is 
in the cone at F z . It now follows from (114) that C = D^, log Hp\{ v i =0 } is diagonal (with respect to the r( 
coordinates) with negative eigenvalues. 

(iv) To complete the proof of the lemma, we now show that the other blocks of the normal Hessian D^^c 
vanish. We firs t no te that since Hp(r z , ■) is constant on Y, the matrix *" vanishes. We also immediately 
conclude from ( 114 ) that 

*' = £> P " A,' logH P \ { ,, r=0} = , 

where p" = (5Rr? m _ r+ i, . . . , 3?ry m ). 



It remains to show that * = 0. Assuming ( |113| ), we expand (110) as before: 

E 7eQ e- (AW ' r> |c A(7) x 7 W| 2 A(7) 



E 



7 eQ l°A( 7 ; 

7n— r 

W0+EMW) + O(hT) 



Hence 



□ 



3.2.3. Proof of Lemma 3_S. It follows from (114) that (Q,w z ) is a c ritical point of ^c{t,',') and hence 
{0} xC z is contained in the critical set. We conclude from Lemma [3.10 that the real rxr matrix B t A~ l B = 
(iB) t (—A~ 1 )(iB) is positive definite, and hence 



det^N^c = det A de^A^B) dctC f , 



(115) 



if z is not a transition point. Therefore {0} x C z is a nondegenerate critical manifold of whenever z is 
not a transition point. 
Recalling (|97|), we set 



b(z) 



-* c (r z ,0,w 2 



-plogH^(-T z ,z) - logH P (T Zj w z ) 



(116) 



We now show that the maximum of SR^c on the r-contour is attained on {0} x C z and hence the maximum 
equals — b(z). We recall that by (jS7|), it suffices to show that SR^t, 0,w) attains its maximum at w = w z , 
and by (|97|), this is equivalent to establishing that log H(t, •), whose Hessian is C, takes its maximum at w z . 
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We first suppose that z is not a transition point. Note that by the above, H(t, •) is constant and has a local 
maximum along Y — fj, F (F z ). Suppose on the contrary that there is a point w € Mp such that 



H(t,w) — sup H(t,w) > H(t,w z ) 

weMp 



(117) 



Let F denote the face that contains (j,p(w). Since H(t, •) is constant on F z , by (117) the face F is not 
contained in F z , and hence C F (~l C° F = 0, where C F denotes the nonboundary points of the cone Cp. Since 
— t G C° F ^ (by the assumption that z is not a transition point), it follows that — r ^ C p . Then we conclude 
from (114) that at least one of the eigenvalues of C(w) is positive, contradicting (117). On the other hand, if z 
is a transition point, then — r G C F ,, where F' is a face whose closure contains F z . Then H(r, w z ) = H(t, w') 
for w' G .F' and H(r, •) attains a local maximum along F'. Then by the above argument, H(t, w') is a global 
maximum. Recalling (116), we then have 



-b(z) = ^ c {t Zi 0,w z ) = sup $c(t z ,0,w z )= sup ^c(r z ,f, w z ) . 
weMp ipeT m ,weMp 



This completes the proof of Lemma pM 



(118) 
□ 



3.2.4. Asymptotic expansions on the classically forbidden region. We continue the proof of Proposition 3.1 (ii). 
Equation (118) says that '5 + b(z) is a phase function of positive type that takes the value on the critical 
set {0} x C z . If z is not a transition point, the points of {0} x C z are nondegenerate critical points. As 
before, we obtain by the method of stationary phase an asymptotic expansion 



K(z) = N^e- Nb ^[ C F(z) + cf (z)N- 1 + 4{z)N~ 2 + ■■■} 



(119) 



valid over each open region 1Z° F . (Since the critical submanifold {0} x C z has codimension 3m — r, the leading 
term of the expansion of e Nb< - z ^ K(z) contains ]\i 2m ~k( 3m - r ) = N 111 ^ .) To show that Hi^p(z, z) also has the 
asymptotic expansion (119), we must show that 



S N {z) = 0(N- k e- Nb ^) V k G N , 



(120) 



where the bound is locally uniform in C 3 for all j. We deform the T m integral in ([71]) to an integral over 
{|CI = T z} as before, so that we have 



S N (z) 



(2^1 [ [ £ e<Q,Ti 

(^) JMp JMp Jt™ \\a\<p 



N 



N 



e -{a,T z +iip) 



fh a (r))m^(w) dipdwdij, 



(121) 



where Rn 



(Np+r. 



(Np)\ 



-i?Ar is rapidly decaying. We first note that 



J2 e^'+^lm^iz 

\a\<p 



< J2 e {a ^\mf(z)\ 2 =H pS (-r z 

\a\<p 



By the Cauchy-Schwartz inequality and the fact established above that H p (t z ,w) takes its maximum at 
w = w z , we have 



a£P 



< J2e- {a ^ /2) \fh P a (v)\e-^^\m^w)\ 

ctSP 



\a£P / \a£P 

= Hp(T Zl n)?H P (T z ,w)? < H p (t z ,w z ) . 
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Therefore, we can write 



S N (z) 



(27T) 



M P JM P JT 



R N (w,r ) )e N ^^' w ' ri;z) d^dwd'n , 



where the phase satisfies 



m < log [H p v(-t z ,z)H p (t z ,w z )] = e 



The estimate (12C ) now follows from the fact that Rn is rapidly decaying. 

Thus, we have verified the asymptotic formula (ii) of Proposition 3.1, together with (a) and (b), and we 
have shown that b{z) is given by ( |116| ). 

Remark: It follows from (|72]), (|S7|), and ( |118[ ) that in a neighborhood U CC C* m of a transition point z°, we 
have uniform bounds of the form 

n|jvp(z, z) < CuN 2m e~ Nb ^ . (122) 
But a sharper inequality should hold; for example, we know that if z° £ dAp, then b(z°) = and 

IWz°,z°) < n% r ;(z°,z°) = [p m + o(l)]N m . 
It is an open question whether we have asymptotic expansions of the form (|119|) at transition points. See 



5.3 for a discussion of the issues involved in studying the asymptotics at transition points. 



3.2.5. The decay function. In this section we verify (c)-(e) of Propositio n |Q| (ii) and then prove Proposi- 
tion We first note that when z S Ap, we have t z = z and hence by (116), b — on Ap. Furthermore, 
since z i— > t z is easily seen to be continuous on C* m and C°° on each TZp, the same holds for b(z). 

We now show that b is C 1 by computing its derivative. Recalling that b(z) — ~^c( T z, 0, w z ; z), we have 

~D z b = D T ty c \(T z ,o,w z ,z) ■ D z t z + D w V c \(T z fi,w z ,z) ■ D z w z + D z ^ c \( Tx ,o,w z ,z) ■ 

Here, D z t z and D z w z are only piecewise continuous, being discontinuous at transition points. By (|8g|)-(|9l|) 
and the fact that (0, w z ) is a critical point of ^c( T z, ■, ■), we have 



\(t z ,0,w z ,z) 



-iD v ^c\(r z fl,w z ,z) = 



The critical point condition also says that D w 1 $c\(t z ,o,w z ,z) = 0- Recalling (|97|), we thus have 

D z b = -D z y c (T,0,w;z)\ T==Tz , w=Wz =- P D z logHv(-T,z)\ T=Tz G C°(C™) . (123) 

Therefore, b S C 1 (C* m ), verifying (e). 

To verify (c), we need to find a more explicit formula for the derivative, using the log coordinates Q = 
Pj + i6j = log Zj. We consider the map £ = . . • , £ m ) : C* m — * C* m given by 

£( z ) = e ^/ 2 . z = e r./2+p+W ; ( 124 ) 



so that ^ s o £(z) = ±(i P (w z ). By ( |L23[ ), 



D p b 



-pD p logHv(-T, e")\ T=Tz 
-pD p \og (l + J2e T,+2p3 



We have 



+ pD p log (l + J2 e2 " 3 
2 e ri+2 Pl 2e Tm+2pm 



(125) 



l + Ee^+ 2 «'"'' 1 + E 



2|6 



2|£* 



2^s ° £ 



(126) 



and similarly, D p log (l + e2pj ) = 2 . Hence 

£>„b = 2p(/i E - Ms ° • (127) 

Let z S C*™ 1 \ yip. If z is not a transition point, then the asymptotic expansion holds at z, and as we 
observed at the beginning of §3.2, TLwp(z, z) is rapidly decaying; hence b(z) > 0. Now suppose on the 
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contrary that z is a transition point and b(z) = 0. Then z would be a critical point for 6, which is impossible 
since Db\ z ^ whenever z £ Ap by (127). This completes the proof of (c). 

It remains to verify formula (|l0|). First we establish ([ll]): Fix a poi nt z in the classically forbidden region, 
and make the change of variables pj = aj /2 + log \ zj\. Then by ( 127 ) , 



b(e T * /2 ■ z)) - b(z) = / ~ D a b{e a l 2 ■ z) ■ da = p 



(e a/2 ■ «) - A*e o e(e CT/2 • z)] ■ da 



(128) 



By definition, q(z) = (ip(w z ) = p^o^(z), and hence PHso{;(e' 7 / 2 -z) = q(e a ' 2 -z). Furthermore b{e T ^ 2 -z)) = 
since e Tz / 2 • z = £(z) e cL4p, and formula (pd|) follows from fll2§|). 



Note that the integral in ( |ll[ ) is independent of the path from to r 2 . Choosing the path <r(r) = rr Zl 
< r < 1, we have 

g(e ff / 2 • z) ■ da = [ * q(z) ■ da = (q(z),T z ) . (129) 
Jo 

On the other hand, since D a log(l + \\e a ^ 2 ■ z\\ 2 ) — /is(e CT / 2 • z), we have 

fi s (e^ 2 -z)-da = log(l + Me—/ 2 • z\\ 2 ) - log(l + ||z|| 2 ) . (130) 

Substituting (12S)— (13C) in (|TT|) , we obtain formula (|l0|). This completes the proof of part (ii) of Proposi- 
tion O. 



We note that b fails to be C 2 at transition points. This follows from (127), (124), and the fact that z i— > r z 
is piecewise C 1 with discontinuous derivative at the interfaces between the regions TZp. 



Remark: An alternate derivation of ( |l0| ) is as follows: By an argument similar to the above, we obtain 



logH P (T,w) = / /ip(e cr/2 • w) ■ da , 
Jo 



and therefore 

log H p (t z ,w z ) = (-Hp(w z ),t z ) 
Formula (0) is an immediate consequence of ( |116[ ) and ( |132| ). 
We now prove Proposition |3.2| : By Proposition [OJ(ii) and ( |116| ), 
1 



(131) 
(132) 



TV 



logII| A rp(z, z) — > -b(z) = plogffs(-T z , z) + logH P (T z ,w z ) 



for all non-transition points z. Let 

u N (z) = ±logIL lNP (z,z)+plog(l + \\z\\ 2 ) = ±log (?)e 2<Q ' P> , *e 



Thus for all non-transition points z, we have 



aeNP 



u N (z) -> Uoo (z) := plog(l + ||«|| a ) - 6(«) = plog(l + \\e T ^ 2 ■ zf) + \ogH P (T z ,w z ) , 



(133) 



(134) 



(135) 



We must show that convergence of (135) also holds at the transition points and is uniform on compact 
subsets of C* m . We again use the log coordinates Q = pj + i9j = log Zj, so that 



D p u N --D p \og ^ ( Je< >--— (JVp) 2<a , p) -]v 



MVP , 



(136) 



and therefore ||Z? p mjv|| < 2p. Since {ttjv(z )} converges for any non-transition point z°, it follows that {itjv} 
is uniformly bounded and uniformly equicontinuous on compact sets. Therefore it converges uniformly on 
compact sets in C* m . □ 
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3.2.6. Precise asymptotics on the classically allowed region. In this section, we prove part (i) of Proposi- 
tion EO. Since 



nAT pS (z, z) 

we have by M) and <M), 



dimH (C¥ m ) O(Np)) _ (Np + m 
Vol(CP m ) ~ V m 



ml = Y\_( N P + j) 

3=1 



U\ NP (z, z) = Y[(Np + j) + R N (z) 

3=1 



where 



Rn(z) 



E 

aGN P S\NP 



\Xa\ 



|XaW| S 



l[(N P + j) 

3=1 



"' 

a£Np£\NP 



NpY, 



(*)| S 



(137) 



Since #(iVpE \ iVP) = 0(N m ), Proposition 3.1 (i) is an immediate consequence of (137) and the following 
uniform asymptotic decay estimate for monomials. 

Lemma 3.11. Suppose that K C C* m is compact, and let U C E be an open neighborhood of ^(K). Then 
there exist positive constants Ci — Ci.k.u, A = U such that 



\C l (K) 



<Cie 



-XN 



a 



whenever — G £ \ J7 , AT > 1 , 1 > 



Proof. We first prove the C° estimate. Choose a positive integer fc > 2m 2 /dist (£ \ U, fj,s(K)). We consider 
the finite collection of polytopes 

:= 2mE + (3 C fcE , /3 G Z m ("I (A — 2m)E . 

For each (3 as above, we let bp(z) be the function on C* m given by formula ([Toj ) with P replaced by Qp and 
with p replaced by k. We claim that the C° estimate holds with 



A = 



inf <{ bp(z) :zeK, (3 G Z m n (fc - 2m)£, -j-Q,g cS\(JbO 



Suppose on the contrary that the estimate does not hold with this value of A. Then we can choose an 
increasing sequence of integers {Nj} and a sequence a? G N(T, \ U) D Z m such that 

2 



e XN > sup 



. JV,-S 



(138) 



By passing to a subsequence, we can assume that ft -> a; G E \ U. Choose /3 G Z m n (k - 2m) E such that 



x G interior^ 



and hence ha? G NjQpfor j > 1. Since dist(x,X) > > diam^Q^), it follows that {\Qp) fl if 
Hence by Proposition 3.2 applied to the polytope Qp, there is a positive constant A such that 



Il\ NQ0 (z,z) < Aexp(-^bp(z)N 
We shall use the multinomial inequality: 



< 



< Ae 



fkN 
\ka 



-XkN 



Vz G K, V./V > 1 



(139) 



(140) 



One way to verify (|l4C| ) is to regard (^) as the partition function Pa? (a) for the polytope E; i.e., the 
number of ways to write a = (3 1 + ■ ■ ■ + j3 N , where (3 G E (see §2.1). Now, given ka G ftiVE we first write 
ka — a + ■ ■ ■ + a (k times) and then decompose each a as above, resulting in f ) different decompositions 
of ka as a sum of kN elements of E. This gives us the lower bound for 



fkN\ 

, ka J 



in A140 
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Wc thus have 



, ^ ^ ->/ ( N j 

a? 



\Xa,i\ 



l 2k 



I 2 = l-^ S | 2 



Since ka 3 € NjQp, it then follows from (139) that 



for all z G iiT, which contradicts ( |138| ). Hence the C° estimate holds. 
We now verify the C 1 estimate: We have 



\fh a (z)\ 



N\ - 



h(z) 



2 



1+ z 



Differentiating, we obtain 



D p \m a (z)\ 



N 



Hz) 



2 z 



N 

a + —D p log/i 



|m a (z)| 



TV 

a + —D p log/i 



and hence the C 1 estimate follows from the C° estimate. Differentiating repeatedly, we obtain all the C l 



estimates. 



□ 



The proof of Proposition 3.1 is now complete 



4. Distribution of Zeros 



In order to deduce Theorem |3j from Proposition p.l\ we need to relate E|tvp(Z s ) to the Szego kernel II|jvp. 
The idea is essentially the same as in [3Z1, Prop. 3.1], but neither the statement nor the proof there cover 
the application we need. Hence we provide a more general statement for the zeros of general linear systems 



on compact Kahler manifolds (Proposition 4.1) 



4.1. Expected distribution of zeros and Szego kernels. As in §1.2, wc let L M be a holomorphic 
line bundle over a compact Kahler manifold, and we consider a linear subspace S C H°(M, L N ), which we 
endow with a Hermitian inner product. (The set of zero divisors {Z s : s G S} = ¥S is called a linear system 
on M.) We let be a local holomorphic frame over a trivializing chart U. We choose an orthonormal basis 
{Sj : 1 < j < k} for S and we write Sj = fj&i, over U. Any section s G S may then be written as 

k 

s = {c,F)ef N , where F = (A, . . . ,f k ) , r.F ^c,/ ; . 

i=i 

The Kodaira map is given by: 

0> 5 : M — > CP*" 1 , * s {z) = [fi(z), . . ., f k {z)] . 



(141) 



(We use the symbol in (141) to indicate that the map $5 is a rational map; it is holomorphic if and 



only if S has no basepoints other than fixed components.) Under a change of orthonormal basis, one gets a 
unitarily equivalent map. The 
the Poincar'e-Lelong formula: 



unitarily equivalent map. The current of integration over the zeros of s = (c, F)ef N is then given locally by 



Zs = 



dd\og\{c,F)\ 



It is of course independent of the choice of local frame ex and basis {Sj}. 

We write a section locally as s = /e^ and we write its Hermitian norm as ||s(z)||/t = a(z)~ 
a(z) = ||ei(z)||^ 2 . The curvature form Ci(L, h) of L is given locally by 



(142) 



' \f(z)\ where 



ci(i, h) = 



2tt 



-dd log a 
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Proposition 4.1. Let (L,h) be a Hermitian line bundle on a compact complex manifold M. Let S be a 
subspace of H [M, L) of dimension > 2. We give S an inner product and we let 7 be the associated Gaussian 
probability measure. Then the expected zero current of a random section s £ S is given by 

E 7 (Z S ) = ^-dd\ogU s {z,z) + c l {L,h) 
where Dq is the fixed component of the linear system VS. 

Proof. The Gaussian measure 7 is the standard Gaussian for an inner product on S. Let {Sj} be an 
orthonormal basis for this inner product. As above, we choose a local nonvanishing section of L over 
U C M, and we write 

k 

s = ^c j S j = {c,F)e L , 
3=1 



where = fjBL, F= ...,/&). As in the proof of [3Z1], Proposition 3.1, we then write F(x) — \F(x)\u(x) 
so that |u| = 1 and 

log\(c,F)\ = log|F|+log|(c,u>| . (143) 



By fll42D, we have 



(E 7 (Z.),v>) = ^J. I (9aiog|(c,F)U)d 7 (c) 



77 
TT 

-\m— l,m— 1 / 



[ (ddlog\F\,<p)drf(c) + ^( [ ddlog\(c,u)\drf(c),<p) 
c k n \Jc k J 



for all test forms tp € T> rn 1,m 1 (U). Upon integration in c, the second term of ( |l43|) becomes constant in z 
and the derivatives kill it. The first term is independent of c so we may remove the Gaussian integral. Thus 

E 1 (Z s ) = ^dd\og\F\ 2 = ^dd(logY / \\S J \\ 2 +loga) . (144) 

Recalling that ILs(,z, z) = ll'Sj'O 2 )!! 2 an d that c±(L, h) = :: ^^-dd\oga, the first identity of the proposition 
follows. 

Before completing the proof of the second identity, we recall the definition of the pull-back of a form by 

9 — 

a rational map X — - > Y: we let T be the desingularization of the graph r C X x Y of g so that we have 
g o 7Ti = 7T2, where 7Ti, 7T2 are the projections from r to X, Y respectively. Then for any smooth form 77 on Y 
we define the pull-back current g*n := Tri^iT^r], which is easily seen to have C 1 coefficients with singularities 
on the indeterminacy locus L g of g. 

We ca n sup pose that U is chosen so that there is a tp £ 0{U) with Div(i/>) = D on U. Write fj = ipgj. 



Then by (144) we have 

E 7 (Z S ) = ^dd\og\F\ 2 = ^-^dd\og\ip\ 2 + ^dd\ogJ2 | 5 . 



^IdB\og\F\ 2 = ^Jdd\og\^\ 2 + ^dBlogj: \ 9j \ 2 



s ■ 



The first term equals Dq. By definition, the second term equals <i>^cj FS outside the indeterminacy locus 
Since I$ s has real codimension > 4 and logj^ \gj\ 2 is (pluri-)subharmonic, the coefficients of j^dd^°sJ2 \9j\ 2 
are measures that do not give mass to i.e., the coefficients are in C l (U). Since the same is true for the 
coefficients of $Jw FS , the current equality extends to all of U. □ 

We note that the expected zero current E 7 (Z S ) is a smooth form outside the base point set of S. We shall 
use the following result on simultaneous expected zeros (away from base points). 

Proposition 4.2. Let M be a compact complex manifold, and let (Li, hi), . . . , (L&, hk) be Hermitian line 
bundles on M (1 < k < dimMj. Suppose we are given subspaces Sj C H (M, Lj) with inner products 
and let jj denote the associated Gaussian probability measures on the Sj (1 < j < k). Let U be an open 
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subset of M on which Sj has no base points for all j . Then the expected simultaneous-zero current of random 
sections si £ Si, . . . , Sk <E Sk is given over U by 



E 



fe fc 

<...X7( t 0^si 1 ...,sJ = f\ E 7j(^j) = A 



-1 



2tt 



dd\ogIl s Az,z) + a(Lj, hj) 



Proof. We first note that the expect current is well defined, since for all choices of the Sj the total mass 



of the zero current Z* 



equals the Chern number c\(L\) ■ ■ •Ci(ife), and hence for each test form (p S 



V m-k,m-k(jj^ the function (Z ait ... tah ,<p) is in x ••• x S k ). 

Next we note the following functorial property of the expected zero current: Suppose that p : Y — ► M 
is a holomorphic mapping of complex manifolds and let p*jj be the pulled-back Gaussian measure on 
p*Sj C H°(Y lP *L). Then 



E,* 7j = p*E 7j (Z Sj ) on p-^f/) S p*S j ) 



(145) 



Indeed, ( |145| ) follows immediately from Proposition 4.1 and its proof. As a special case, suppose that Y is 
a submanifold of M and p is the inclusion. Then ( |l45| ) becomes 



E 71 ([y]AZ flj ) = [y]AE 7 ,(Z Sj .) on £/ 



(146) 



(In fact, (146) holds for any subvariety Y, but we don't need this fact here.) 



We now verify the current identity by induction on fc. For fc = 1, this is Proposition 4.1, so assume that 



fc > 1 and the proposition has been verified for fc — 1. Consider Y 



J Sl,...,Sfc_l I 



By Bertini's Theorem, we 



know that Y is smooth (and of codimension fc — 1) for almost all s±, . . . , Sk-i- Therefore by (146), 

z si,...,s k djk(sk) = J Z s ^..., Sk - x A Z Sk dy k (s k ) = Z 8 ^...,s h - X A E 7fc (Z s J for a.a. s X) . .., s fe _i 
Integrating over s±, . . . , s^i and using the inductive hypothesis, we obtain 

E 7ix---x 7fc (Z Sl ,... iSfc ) = E 7lX -x7fc_i(2'ai,...,sh_i) A E 7fc (Z Sfc ) = E 7l (Z Sl ) A • • • A E lk (Z Sk ) . 



The other equalities follow from Proposition 4.1 



□ 



If we consider the case where M = CP" 1 , we obtain a formula for the expected simultaneous zero current 
for random polynomials with given Newton polytopes: 

COROLLARY 4.3. The expected zero current ofk independent random polynomials fj 6 H°(CF m , 0(jpj), Pj), 
1 < j 5= k> *s given over C* m by 



E 



° h ( Z fu-, 



fe 

A 



2tt 



-ddlogllip (z, z) +PjU. 



c?" 



Proof. Since 77°(CP m , O(pj), Pj) has no base points in C* m , the conclusion follows from Proposition 4.2 and 
the fact that c\{0{pj) 1 h FS ) =pjLu FS . □ 



Remark: By Proposition 4.1, the identity of currents 

— 1 



E|p(Z / ) = ^— 3dlogII|p(z,z) 



holds on all of CP m . 

We say that p is minimal for P if 



min ctj = (1 < j < n) , 



max a 

a6P 



(147) 



(148) 



(Otherwise, we can choose p to be smaller, after translating the polytope if necessary.) Th e hyp othesis that 
p is minimal for P guarantees that H° (C¥ m , O (p) , P) has no fixed divisors. In this case, ( 147 ) holds as an 
identity of (1, l)-forms with C\ oc coefficients. 
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4.2. Asymptotics of zeros. In this section, we prove Theorems and and give some additional asymp- 
totic formulas for zeros. Theorem |l| is a special case of Corollaries]! and g of T heorem [|. 

To prove Theorem ^, we let un, u<x> be as in the proof of Proposition 3^. By Corollary h3 (recalling 



that u>„ 



l -dd\og{l + ||z|| 2 )), we have 



-^-ddu N = —E\ NP (Zf) 



and hence un is plurisubharmonic. By (135) and (1491), we conclude that 



— E\ NP (Z f ) 



-1 



2tt 



ddu r 



-1 



2- 



996 



(149) 



(150) 



where differentiation is in the distribution sense. 

We now show that the current ddb 6 7-?' 1 ' 1 (C* m ) is given by a (1, l)-form with piecewise smooth coeffi- 
cients; in fact ddb is C°° on each of the regions TZp given by (||). (Equivalently, the Radon measure ddbAu> FS 
does not charge the set of transition points, and hence formula ( |150| ) can be interpreted as differentiation in 
the ordinary sense on the regions IZp.) By Proposition [O], b e C 1 (C* m ); i.e., if z° e 97?. p n dTZpi, then the 
values of db(z°) computed in the two regions 7?.i? and 1Zf> agree. Then for a test form ip G 2? m_1,m_1 (C* m ), 
we have 



(596, <P) = V. [ bddip = V / ddb/\<p-J2 [ 



(96 A + 6 A 9^) 



(151) 



Note that 97?. f consists of C°° real hypersurfaces (consisting of those points of dlZ p that are contained in 
the boundary of only one other region TZ pi) together with submanifolds of real codimension > 2, and hence 
Stokes' Theorem applies (see e.g. 4.2.14]). Since 6 and 96 are continuous on C* m , the boundary integral 
terms in (151) cancel out and we obtain 



(996, = 996 A ip = f 

p Jn% Jc 



ddb A ip 



(152) 



C*™\E 



where E := \J F &R.p is the set of transition points. Formula (152) says that the current 996 is a (1, l)-form 
with piecewise smooth coefficients obtained by differentiating 6 on the regions 1Z° F . 
We now let 



2tt 



-ddu c 



puj ¥ 



2tt 



996 



(153) 



on each of the regions TZp. By ( 152 ), the current ipp is also a piecewise smooth (1, l)-form; by ( |15C 

7V~ 1 E| A r P (Z / ) -> %p P (weakly) . (154) 

We shall show C\ oc convergence of (154) when we prove Theorem || below. Continuing with the proof of 
Theorem ||, we observe that (ii) is an immediate consequence of (153), since 6 = on the classically allowed 
region. 

To verify (iii) , we again use the log coordinates pj + iOj = log Zj , so that 

Woo =piog (i + y]i 

Since depends only on (pi, . . . , p m ), we have 

~ T 



b(e") 



(155) 



2tt 



: 99u„o = 



dQ A dC,k > 



(156) 



Thus we must show that the Hessian of u oa (p) has rank r at points z° € TZp, where r = dimF. From 
(^25|)-(p6l) and ( ||| ), we obtain 



2p Ms o £ 



(157) 
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where £(z) = e T */ 2 • z as before. Hence ^PpUoo equals the Jacobian of /is ° £■ Since /is ° £,{ z ) = h^p(w z ), 
we easily see that 



(158) 



so rankD(/is°£) < r. In fact, (158) is a submersion, so that the rank equals r. To see that it is a submersion, 
we recall from the proof of Lemma 3.7 that the 'lipeomorphism' <I> -1 : S° — > 7V° restricts to a diffcomorphism 

S° D Me(^f) " F x C F C N° , ^ (|mpK), - (/is o £(z), -t z ) , 

and hence (158) is a submersion, completing the proof of (iii). 

We now prove Theorem f|; the case fc — 1 of the theorem will then yield part (i) of Theorem ||. 
Let Pi, . . . , Pfe be Dclzant polytopes, as in Theorem |[ We first show that 



N-% NP {Z fu ... Jk ) — * f/'fe := ipPi A ■ ■ ■ Atpj 



(weakly) . 



(159) 



For 1 < j < k, we let 



,3 - 



W N = logU\ NPj (z,z) +pjlog(l + ||z|| 2 ) , 



so that, recalling ( |135| ), 

v? N {z) - <,(z) := ft log(l + ||z|| 2 ) - = ft - log(l + \\e T f 2 ■ z\{ 



log H Pj {t 3 z ,wD 



where t° z , w 3 z are as in Lemma with P — Pj. By ( ]153| ) , 



2tt 



Recalling Corollary 4.3, we introduce the (k, fc)-forms 



-1 



2tt 



A • • • A 



-1 



2tt 



99< I A ; ' K v/ , Vft (Z; /J 



Since 



'AT 



it jo locally uniformly, it follows from the Bedford- Taylor Theorem |BT, KJ that 



K N 



-1 



2tt 



<9<9u:L ) A 



2tt 



(weakly). 



(160) 



(161) 



In fact, the limit current in (161) is absolutely continuous, and hence is equal to the locally bounded (piecewise 
smooth) (k, fc)-form ipk- To see this, we note that by the Bedford- Taylor Theorem, 



k 

A 



k 

A 



where tp e denotes an approximate identity. Since the currents ip Pj have coefficients in £j^ c , the forms 
Aj=i(V'p * Pe) have locally uniformly bounded coefficients; absolute continuity of the limit current follows. 
The weak limit (159) now follows from (160)— (161). 

To complete the proof of Theorem f|, we must show that kn — > ipk in ^(K) for au compact K C C* m . 
Let e > be arbitrary, and choose a nonnegative function rj S 2?(C* m ) such that 77 = 1 on a neighborhood 
of (Pi U • • • U Ek) H if, where Pj is the set of transition points for Pj, and 

ipk A r]uj r : l - k < e . 



By fll59D-(ll6C|), 

/ hn A ?/w™" fe < 2e for TV > . 
(Note that the above integrands are nonnegative.) Therefore, 

\\*N ~ 1pk\\ci(K) < IK 1 - V)(*N - i>k)\\ c i(K) + \\VKn\\ci( K ) + \\vi>k\\ c i(K) ■ 
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Since kn is a positive (k, fc)-form, 

Similarly, \\vi>k\\ c i(K) = Ik ^ A 7 ? w "s~ fc < e - Since 

(1 — ?j)(«;jv — ^fc) — * uniformly on if , 

it follows that 

\\k n - ipkWc^K) -> - 

completing the proof of Theorems || and [|. □ 

We note that the conclusion of Theorem |J can be replaced with an asymptotic expansion away from 
transition points: 

Theorem 4.4. Let P\, . . . , P^ be Delzant polytopes. Let U be a relatively compact domain in C* m such that 
U does not contain transition points for any of the Pj . Then we have a complete asymptotic expansion of 
the form 

j^\NP u ...,NP h {Zf 1> ...j h )^^p 1 h---hiJp h + ^ + --- + j^ + --- on U, 
with uniform C°° remainder estimates, where the (pj are smooth (k, k)-forms on U. 



Proof. By (|134|) and ( J160| ) we have 

1 k V—l 

w E lNPu ... tNPk (z fu ... Jh ) = f\ —dd 



— logIL lNP .(z,z)+pjlog(l + \\z\\ 2 ) 



The conclusion now follows from the asymptotic expansion of Proposition 3.1. □ 

The proof of Theorem |3|(iii) gives us some more information about the expected zero current in the 
classically forbidden region, which we state in the theorem below. We first define the complexified normal 
cones 

C F = {T + i9 : T e C F , 6 S Tp-} . 

(Recall that Cf C Tp.) We note that Cf is a semi-group, which acts on TZp by the rule j](z) — e v ■ z; 
we call this action the '(joint) normal flow.' The (maximal) orbits of the normal flow are of the form 
Cf ■ z° = {e v ■ z° : n £ Cf}, where z° £ ^^(F). We note that the orbit Cf • z° is a complex (to — r)- 
dimensional submanifold (with boundary) of C*"\ (Indeed, (Cf- z°)nlZ° F is a submanifold without boundary 
in Tip.) 

THEOREM 4.5. Let P be a Delzant polytope and let tpp be the limit expected zero current of Theorem |^. 
Then ipp vanishes along the orbits of the normal flow. 

Proof. Let 

O := C F ■ z° = {e T+l6 ■ z° : r € C F , 9 € Tp} 
be a maximal orbit of the normal flow, where /is(z°) £ F. For z — e T+t6 ■ z° £ O, we have 

A^s ° £(z) = -p-p{w z ) = -fi P (w e r. z o) = -fi P (w z a) = nz(z°) 

p p p 



and hence /i^ o £ is constant on O. It then follows from (156)-(157) that iJjp\o = 0. □ 



Remark: Theorem 4.5 tells us about the behavior of tangent directions to typical zero sets \Zf \ of polynomials 
/ £ H (CP m ,O(Np),NP) as N — > oo. Roughly speaking, the tangent spaces of \Zf\ are highly likely to 
be close to containing the tangent spaces of the orbits of the normal flow, for large N. E.g., let us consider 
a region 7Z° F , where F is an edge (i.e., dimF — 1). Let T^ ow denote the holomorphic tangent space to 
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the normal flow through a point z € 1Z° F . Then dimcTj iow = dime Zf = m — 1, and for any compact set 
K c 1Z F , we have 



^E| WP ^ ^dist(T^° 2 , T*° w ) 2 <*Voi 2/ (z)^ 



where 'dist' means the distance in the projective space of complex (m — l)-dimensional subspaces of T ( 



(162) 



To verify ( 162 ), we let V be a compactly supported vector field on 1Z° F such that for all z £ 1Z° F , V z is a 
(1, 0)-vector tangent to the normal flow through z. Let ip = *(iV AFje P m ~ 1,m_1 (7?.^) be given by 



AA^ = (A,^A7 2 )> p m 3 



By Theorem 4.5, we have 



^p, <p)= Nj p Mp= / (fa, iV z A iw™ = 



we let r f ec°°(z]; cg ,r™ 



x ; y " ~ r 1 1 1 ) denote the dual to the volume form on Z T f s , so that 



\\T f A V A V\\ dVol Zf 



(T f ,v)dVol z , = (Z f ,<p) 



Hence by Theorem |](i), 
1 

NT 

which yields (|l62| ). 



E 



\NP 



\\T f AVA V\\ <No\ Z} = -E lNP (Z f , <p) -> (i/j p , <p) = 



(163) 



(164) 



As mentioned in the introduction, we can apply Proposition |4.2| (as in [SZ1 1 for the case k = 1) to the 
asymptotic expansion of the Szego kernel in pc| to immediately obtain the following asymptotics using toric 
norms: 

Proposition 4.6. Let (Mp,Lp) be as above, and let Ij,6j be the action-angle variables of the moment map 
fip : Mp — > R m of the T m -action on Mp. Then we have: 



-L Ei M P (z fu ... tfk ) 



Proof. Let wp = X^li ^ d6j be the natural Kahler form on Mp. We recall the asymptotic formula for 
the Szego kernel in [Kq]: 

1^»^ FS = uj p + 0(1/N) . 
By Proposition 4.2 with M = Mp, L = Lp, we then have 

E Mp (Z Sl ,... iSk ) = (E Mp (2 s )) fc - = A fc [ w * + 0(1/JV)] . 

□ 

4.3. A moebas in the plane. The term 'amoeba' was introduced by Gelfand, Kapranov and Zelevinsky 
|GKZ | to refer to the image un der the moment map of a zero set Zf lt ...j k of polynomials, and have been 
studi ed in various contexts (see [ FPT , GKZ, Mil , PR| and the references in the survey article by Mikhalkin 
[ Mi2 1 ) . The image of a zero set under the moment map /is is called a compact amoeba, while the image 
under the map 

Log :C™->R m , (z 1 ,...,z m )^(log|z 1 |,...,log|z m |) 
is a noncompact amoeba, or simply an amoeba. Note that Log is the moment map for the T m action with 
respect to the Euclidean sym plectic form ^ dxj A dyj, and Log = hC o fx^, where C : S° w R m is the 
diffeomorphism given by (103). 

To illustrate what our statistical results can say about amoebas, we consider zero sets in C* 2 . An amoeba 
in R 2 is the image of a plane algebraic curve under the Euclidean moment map Log. An example of an 
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Figure 5. An amoeba with polytope 4S 

One notices that this amoeba contains 12 'tentacles'. By definition, a tentacle on a compact amoeba A 
is a connected component of a small neighborhood in A of A n 9E; the tentacles of a noncompact amoeba 
correspond to those of the compact amoeba under the diffeomorphism E° as R m . 

There is a natural injective map from the set of connected components (which are convex sets) of the 
complement of a noncompact amoeba A to the set PC\1? of lattice points of the polytope, and that there are 
amoebas for which each lattice point is assigned to a component of the complement. (This fact is also valid in 



higher dimensions; see [FPT, Mil].) For a generic 2-dimensional noncompact amoeba with polytope P, each 
lattice point in dP corresponds to a distinct unbounded component of M 2 \ A, and adjacent lattice points 



correspond to adjacent unbounded components. (The correspondence is given by [Mil, §3.1] or |FPT|.) 
Hence the number of tentacles of A equals the number of points of dPDZ; this number is called the length of 
dP. (For example, the 12 tentacles in Figure || correspond to the 12 lattice points of <9(4E).) Each tentacle 
corresponds to a segment connecting 2 adjacent lattice points on dP. 

We can decompose dP into two pieces: d°P = dP flpS and d e P = P D d(pE). Tentacles corresponding 
to segments of d°P end (in the compact picture £) at a vertex of E, and tentacles corresponding to segments 
of d e P are free to end anywhere on the face of E containing the segment. We call the latter free tentacles, 
and we say that a free tentacle is a classically allowed tentacle if its end is in the classically allowed region 
Ap. For an amoeba A, we let vat(A) denote the number of classically allowed tentacles of A. It is clear 
from the above that 

vat(A) < #{frce tentacles} = Length(<9 e P) 

and that this bound can be attained for any polytope P. Here, 'Length' means the length in the above 
sense; i.e., the diagonal face of pE is scaled to have length p. As a consequence of Theorem [| (for m = 1), 
we conclude that this maximum is asymptotically the average: 

COROLLARY 4.7. For a Delzant polytope P, we have 

^E| A r P (z, AT (Log(Z / )) ) -> Length(3 e P) . 

Proof. Let Fi,F%, F3 denote the facets of pE and apply Theorem ||(ii) to the 1-dimensional polytopes PHFj, 
3 = 1-2,3. ' □ 



2 The authors would like to thank T. Theobald for giving us permission to use this figure from his paper. 
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5. Examples 



We now compute b(z) and ipp for some examples of toric surfaces (m = 2). Recalling ( |10|) and (135), 

1 + IMP 



b(z) = q(z) ■ t z +plog 



q(z) =p^(e T * /2 -z) G dP 



1 + \\e^/ 2 ■z\\ 2 / 
b(z) +plog(l + ||z|| 2 ) =plog(l + \\e T */ 2 ■ z\\ 2 ) - (q(z),r z ) 



Hence by (|153| ), 



2tt 



-ddu r 



2n 



-dd 



p\og(l + \\e^/ 2 -z\\ 2 )~(q(z),T z ) 



(165) 



5.1. The square. For our first example, we let P be the unit square with vertices {(0, 0), (1, 0), (0, 1), (1, 1) } 
so that p = 2 and M P = CP 1 x CP 1 . Recalling that 



Ht,{zi,z 2 ) = 



Z2 



i + M 2 + N 2 ' i + Uil 2 + U 2 | 2 



we see that the classically allowed region is given by 

|*i| 2 -i<H 2 <N 2 + i, 

as illustrated in Figure [j] in the introduction. The forbidden region consists of two subregions: 
Tl F = {{z u z 2 ) : M 2 > \ Zl \ 2 + 1} , F = {(an, \) : < x x < ±} , 

Tl F . = {(Z!,z 2 ) : \z 2 \ 2 < \ Zl \ 2 - 1} , i* 1 * = {(±,z 2 ) : < x 2 < i} . 

Suppose that z is a point in the upper forbidden region TZp. Write r z = (n, t 2 ); then n = since r _L T^?. 



Let 

Writing 
we have 

Therefore 



o(z) = 2^ E (e^ /2 • z) = (a, 1) G 2F C <9P 



N 2 =«i, N 2 =s 2 , |e r2 / 2 z 2 | 2 =e r2 s 2 = s 2 , 



si 



A - 2 



1 + Si + S 2 2 1 + Si + §2 2 



1 



•Si 



1-a 



s 2 



si si 
a 



1-a a' " 1 + si l + |zi| 2 ' 
si l + ki| 2 



s 2 /s 2 



We have 



(9(2),^ 
Therefore 

We conclude that 



log(l + ||e^/ 2 .z|| 2 )=log(l + |z 1 | 2 + >|L> 

1 + Nil 2 



as 2 \Z2\ 



i + N 



1 , 0,lo. 



^2 



= log(l + |zi| 2 ) + log2 
log(l + |zi| 2 )-log|z 2 | 2 . 



uoo = log|z 2 | 2 + log(l + |zi| 2 ) +log4 



ip P = < 



^aaiog(l + |zi| 2 ) for |zi| 2 + l < |z 2 | 2 

2cu FS = ^dd\og(l + \ Zl \ 2 + \z 2 \ 2 ) for \ Zl \ 2 - 1 < |z 2 | 2 < | 2l | 2 + 1 



I ^d31og(l + |z 2 | 2 ) for \z 2 \ 2 < N 2 -l 

(where the third case is by symmetry). Note that ip P has constant rank 1 in both of the forbidden regions 
IZp, Hf*, as indicated in Theorem H(iii). 
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By the above, we also obtain 



-b(z 



4|z 2 | 2 (l + |zi| 2 ) 



for |zi| 2 + l< \z 2 \ 



(l + | Zl |2 + |z 2 |2)2 

and similarly for | z 2 1 2 < | z\ | 2 — 1 . 

Remark: On the boundary {|zi| 2 = |z 2 | 2 — 1}, we have e~ b ^ = 1 as expected. On {\zi \ = c}, we have the 
growth rate e~ b ^ ~ 1/U 2 | 2 as z 2 — > oo. 



5.2. The Hirzebruch surfaces. We now consider the trapezoidal polytope of Figure f| below. 




Figure 6. P = {(0, 0), (1, 0), (2, 0), (0, 1), (1, 1)} 

Again p = 2, but this time Mp A P 2 is the blow up of (0,0,1), i.e., Mp is the Hirzebruch surface F\ 
(see [JFuJ). Comparing with the case of the square, we see that the classically allowed region is given by 
\ z 2\ 2 < l£i| 2 + l, and the forbidden region coincides with the upper forbidden region TZf = {| z 2| 2 > l z i| 2 + 1} 
from §5.1. Thus, the map z <— > (t z ,w z ) is the same as before when z is in the forbidden region TZf- Hence, 
g-K 2 ) and ijj p are a is the same as in §5.1 on TZf- On the classically allowable region, e^ b ^ = 1 and 
i>p = 2uj fs . 

5.2.1. The Hirzebruch surfaces F n (n>2). Let n > 2 and consider the polytope 

P = {(0,0), (1,0), ...,(n+l,0), (0,1), (1,1)}, 
so that Mp — F n (see jF\J). Here p — n + 1 and ip has an interior vertex v = (^j, ^j) G S°. 
We see that the classically allowed region fJ-^i^P ) is given by 

H 2 + i _J_ 

n n — 1 

This time, the forbidden region consists of three subregions: TZf, TZv, TZf 1 , where 

F = {( Xl , jjij) : < x x < ^} , F' = {(x!,x 2 ) : a; 3 =I(l -an), ^ < a* < 1} 

(See Figure fj] below.) 

Suppose that z is a point in the region TZp. Then t z = (0,t 2 ), as before. Let 

q(z) = {n + l)Ms(e Ti/2 • z) = (a, 1) € (n + 1)P C <9P (0 < a < 1) . 

Again writing 

N a =*i, = s 2 , |e r2/2 z 2 | 2 =e T2 s 2 = s 2 , 

we have 

si a s 2 1 



|z 2 | 2 < min 



§2 71+1 1 + Si + S 2 Jl + 1 
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n + 1 




n + 1 



22 



Figure 7. P = {(0,0), (1,0), . . . , (n + 1,0), (0,1), (1,1)} 



Therefore 



In particular, 



and therefore 



We have 



si 



n — a 



S 2 



1 Sl 



— , a = 



nsi n\Zi\ 



n — a a 1 + S\ 1 + \zi\ 2 ' 

si 1 + Ui| 2 



e r2 = s 2 /s 2 = 



aS2 «P2| 

1 



a <1 ^ \zif < 



R-j- = {(zi,z 2 ):\z 2 f>^- +1 '~ 12 



n- 1 



n-l 



f IK :||") = log ( 1 - hil' + ^SW) = LosCI + |: ^ ) ■ log 



n|z 2 

2 



n+ 1 



+ N 



log(l + |zi| 2 ) - log |z 2 | 2 - logn . 



Therefore 
Hence, 



n\z 2 y 

Uoo = log |z 2 | 2 + nlog(l + |zi| 2 ) + (n + 1) log(n + 1) - nlogrc . 



V>p = n^— ddi og (i + M 2 ) 

zi7T 



-6(z 



2\n 



} (n + l)" +1 |z 2 | 2 (l + | Zl | 2 ) 



> for z G . 



n» (1 + N 2 + N 2 )" +1 J 
Now suppose that z is a point in 7\Lf<. Since r z _L Tp>, we can write r z = (n,nri). Let 



-J-j-g(z) = M S (e W2 ■ 2) = (a*, -(1 - Xl )) € F> . 
n + I \ n J 



As before, we write 



si = |zi| 2 , s 2 = \z 2 \ 2 , 



h = \e T ^ 2 Zl \ 2 = e Tl Sl , S 2 - \e nT ^ 2 z 2 \ 2 = e nT 's 2 . 



(166) 
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By (166), we have 



Solving for s*i, §2, we obtain 



si 



1 + si + s 2 



A'2 



1 



Si 



Xl 



1 + si + §2 n 
1 



n — 1 1 — xi 



S2 



n — 1 



(167) 



Therefore, e nri = S2/S2 = | ^2 1 2 /(n — 1), so we have 



Thus, 



Ti 



log(l + \\e T > /2 ■ zf) = log(l + Si + S 2 ) = log 



log(n-l)|* 2 | 2 



si 



(n-lJVn^al 3 /" 



Ti-l (n- iy-l n \z%\ 2 f r - 



By fll66|), 



(q(z),r z ) = (n + 1) ^f^ 1 ' ~( x ~ ^i)^ ' ( T i' nr i)^ = (« + 1)ti = ~"~~~~" log ( n ~~ ^N 5 
Hence by fll~65l), 



V^p = (n + 1)— — 99 log 

Z7T 



for z G 7£p' 



4 n-l (n- l) 1 /»| 22 |2/« i 

(Note that has constant rank 1 on lZp> as indicated by Theorem [|(iii), since ipp = (n+ l)g*uj CP i on IZpi , 
where <? is the multi-valued holomorphic map to CP 1 given by g(z\, z%) = {c n Z\, z\ ).) 

By Theorem [|(iii), we know that ipp — on 1Z V . To complete the description of ipp, it remains to describe 
the regions TZf' and 1Z V . We note that a forbidden point z lies in TZf' if and only if x\ > jjxj. By (167), 
this is equivalent to si > ^zy, or 

N 2 < (n-l)"- 1 ^! 2 " . 

Therefore 



Kf> = \ (zi,z 2 ) 



' <|z 2 | 2 <(n-l)'^ 1 |z 1 | 2n , \ Zl \ 2 > ' 



This leaves us with 



71 — 1 



K v = ( Zl ,z 2 ) : M 2 > («- l)"" 1 !^! 2 ", Nil' > 



n — 1 



n - 1 



To summarize: 



(n + i)^i 9 aiog(i + |z 1 | 2 + |z2i 2 ) 

rPgc>dlog(l + |* x | 2 ) 



for |z 2 | 2 <min{M!±i j _l T | 
for N 2 >^, ki| 2 <^r 



(n+1) ^gl 0g (^ + __^ll_ ) fa. _i_< N 2 <(n _ ir -l kl| 2„ 

for |z 2 | 2 > (Ji-l)"- 1 !^! 2 ", |zi| 2 >^t 



5.3. Remarks on the behavior at transition points. So far, we have studied the uniform asymptotics 
of TL\ N p(z, z) on compact subsets of the classically allowed region and on compact sets without transition 
points in the forbidden region. We now address a few words to the asymptotics at the boundary of the 
classically allowed region, which is its caustic set. For background on oscillatory integrals and caustics, we 
refer to ]AGV| ]. 

We thus consider Ui NP (z, z) as an oscillatory integral (denoted K N (z) in ([72])) on (C*) m . Since it is T m 
invariant we may regard it more simply as an oscillatory integral in z — e p on R™ or on pE, with the phase 
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W of ( |73| ) regarded as a function on T m x Mp x W+. The asymptotics of H\np(z, z) are determined by the 
component 

C° = {(0,«;; e") € x M P : ^jj(e") = MpH}, (168) 

of the critical set on which tp — 0, so we may (and will) ignore other critical points which may occur. The 
map (0, w; e p ) — ► w gives a natural identification of C° with Mp. 

As is standard in the theory of oscillatory integrals depending on parameters, we define the map 

i* : C*° -> T*M™ 4*(0, w\ e p ) = (e p , d„¥), (169) 

and we denote the image by A = t*(C°). Over the set of parameters z for which the critical point set C° 
is a non-degenerate critical manifold, i-q, is a fibration and A is a Lagrangean manifold. The set of z = e p 
for which the phase is degenerate is known as the caustic set. Under the identification of C° ~ Mp, we have 
id z ^(w; e p ) e on C°, hence A is the zero section of T*P. Identifying with the map [ip x 0, we find 
that the caustic set is the boundary of the classically allowed region. 

Uniform asymptotic expansions of IL\np(z, z) in the vicinity of the caustic could be obtained with a careful 
study of the singularities of the moment map fip. The uniform expansion would give a smooth transition 
in a boundary layer from the N m mass density in the interior of P (the illuminated or allowed region) to 
the e~ NhtyZ ^ mass density in the exterior (the shadow or forbidden region). The order of the asymptotics 
along the caustic depends on two competing effects: on the one hand, oscillatory integrals whose (analytic) 
phase has a degenerate critical point decay more slowly than in the non-degenerate case (by the index of 



the singularity at the critical point, see AGVQ. On the other hand, the principal term of the amplitude 
may vanish on the caustic, causing more rapid decay. We will be content to illustrate this in the simplest 
example. 

5.3.1. Example: P = [0, l],p = 2. We put m = 1, P = [0, l],p = 2,p£ = [0,2]. Thus, Mp = CP 1 . The only 
interior facet is F — {1}; recalling (|73]), our phase is: 

*( w) = 1 o g i±i^ +2 ,„ g i±i«. (170, 
1 + \w\ z 1 + \z\ l 

We are interested in the asymptotics of TLmp(z, z) when 

2|z ' 2 

We claim that IIijvj>(,z, z) ~ CN 1 when \z\ = 1. 
i i 2 

Since fip(w) = .^J i 2 , we have fip(w) = {1} w = oo. To determine the asymptotics of TL\np(z, z) 

at \z\ = 1 , we put r] — ^ S C and consider the Taylor expansion of the phase at 77 = 0, tp = 0. We change 
to polar coordinates i] = re t9 around the fixed point of the toric S 1 action. The phase and amplitude are 
independent of 9 and the principal part of the Taylor series (at \z\ = 1) has the form: 

W 1) ~ i W 2 - \p 2 + 0(l^| 3 + \^\ 2 V? + l^lkl 4 ) • 



^(z) = = 1 



By a general result [ AG V 1 , the asymptotics of (generic) oscillatory integrals are determined by the principal 



part of the Taylor expansion of the phase at the critical point, which gives the model integral 



TV />00 



I N = / e N ^ r -X* ^A(r,(p)rdrd(p, 

J -IT JO 

where A is an amplitude decaying sufficiently rapidly at infinity. Here, we also used that the amplitude 
of H\np(z, z) is simply the volume density in our local coordinates (</?, r, 9). We observe that the phase 
fir, (p) = iipr 2 — j ip 2 is weighted homogeneous in the sense that 

f(N 1 /\N 1 / 2 <p) = Nf(r,tp). 
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We also observe that the amplitude vanishes to order one as r — » 0. We change variables as indicated to 
obtain that 

U\ NP (1, 1) ~ N 2 I N = N 2 [ [ N^e^+c^^piN- 1 / 4 ^ N~^ 2 ip)rdrdip ~ CN , 



as claimed. 

Alternatively, we observe that the index of oscillation of the phase with \z\ = 1 equals —3/4. This agrees 
with the general result (under assump tions s atisfied here) that the index of oscillation equals the 'remoteness' 



of the Newton polygon of the phase (|AGV|, §11.6, Theorem 4), i.e. the reciprocal of the value oft such that 
the ray (t, t), t > intersects the boundary of the Newton polygon. The index of the singularity is therefore 
S = 1/4, and therefore I N ~ CN- 1 . 
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